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UPPER AND LOWER FUNCTIONS FOR 
DIFFUSION PROCESSES 


S. K. ACHARYA AND M. N. MIsHRA 


Department of Statistics, Sambalpur University, Jyoti Vihar, Burla, Sambalpur 
Orissa 768017 


(Received 18 August 1987; after revision 19 January 1988) 


The paper is concerned with the study of upper and lower functions for dif- 
fusion processes described by the non-linear time-dependent stochastic 
differential equation. 

1. INTRODUCTION 


Let X (t) be the solution of the non-linear time dependents tochastic differential 
equation 


dX (t) = a(X (t), t) dt + dW (t) (1.1) 
with the initial condition 
X (0) = X, me ee 
where W (t)is a standard Wiener process and X, is independent of F {W(t), t > 0} 
with E(x ) < co, 
Let A > 0. Let Hy, betheclassof non-negative, non-decreasing functions defined 
on [A, co) which increase to oo with f. 


For h(t) € Ha, we say that / (t) belongs to the upper class or lower class acco- 
rding as 
P{X(t) >h(t)i. 0. as t > c0} = Oor 1. 
In this paper our aim is to develop the integral test criterion for the solution pro- 
cess of (1.1) to decide whether h (t) belongs to the upper class or lower class. 


We give preliminary lemmas in section 2. [n section 3, Theorem 3.1 gives a _ re- 
sult analogous to Strassen’s invariance principle. The integral test criterion for dif- 
fusion processes described by equation (1.1) have been developed in Theorem 3.2. 


Problems of the above type have been considered by many authors. Let {X,} be a 
sequence of independent random variables. Feller? and Chung" have studied the asym- 


ptotic growth rates of S, = >» X,and M,= max =(|5S;| respectively, which are 
i=] 1<ixn 
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considered to be fundamental papers in this area. Incase of Brownian motion {W (t), 
t > 0}, Kolmogorov has developed an integral test for non-decreasing function h (t) 
so that A(t) belongs to upper class or lower class according as the integral converges 
or diverges. The same problem has been considered by Strasson® for a martingale dif- 


ference sequence {Y,} and by Jain et al.® for the partial sum S, = at Y;. Jain and 


Taylor’ have studied the asymptotic growth rate of M '(t) = Pie |W (u) |. For 
t 


the first time Mishra and Acharya® have developed the integral test criterion to decide 
whether A (t) belongs to the upper class or lower class for diffusion processes described 
by the homogeneous stochastic differential equations of the Ito type. 


2. NOTATIONS AND PRELIMINARIES 


This section is devoted to the background materials which have been used in this 
paper. Let us consider the stochastic differential equation (1.1) where W (t)isa 
standard Wiener process. We assume that a (x, 1) is reelvalued, well defined, measur- 
able for x € (— ©, co) and satisfy the following conditions, 


(A,) for some constant K and0 <8 < 1, 
la(x,t)| << K/Q+ |x] yrs 

(A,) force > 0 and x, y in (— ©0, ©0), there exists a constant Le such that 
|a(x,t)—a(y,t)| €Le|x—y| 

where | x! <cand|y|<ce. 


Under conditions (A;) and (A,) and the initial condition XY (0) = X,, Gikhman 


and Skorohod‘ have shown that, there exists a unique solution X (t) of (1.1) in an arbi- 
trary time interval [0, 7] and 


X(th= X, + a(X (s), 5s) ds +W(t). Be BB 


Moreover X (t) is a Markov process whose transition probability is given by 
P (Xo, t, A) = Px, (X(t) € A), eZ) 
To prove the main theorem we need the following lemmas. 


Lemma |—Under condition (A,) 
t 
E |J a(X(s), s)ds}* = O (21-8), 
0 


For the proof of this Lemma refer to Friedman? (p. 184), 
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Lemma 2—Let ¢ (t) increase monotonically to infinity with t and {W1t), t > 0} 
be a Brownian motion process. Then 


P{W (t) > UP ¢(t)i. 0. as t + cc} = Oor 1. 
according as 


1 


is convergent or divergent. 
The above result is due to Kolmogorov [see Ito and McKean’, p. 165]. 


Lemma 3—Let g be an eventually non-increasing function from [0, ce) to [0, eo) 
and | be a measurable function from [A, ee) to [0, co), forsome fixed A > 0. For 
h € Aa, define 


F (h) = Ps (h() v(t) at 


which may be either finite or infinite. Assume that 


(a,) for every h € Ha and for every B such that 
B 

B>A>0O, f g(h(t)) ¥ (t) dt < 9. 
A 


(a,) There exists 4,, h,, two members of Ha, such that 


h, & h., F (hz) < 0°, while F (41) = 9 and 
B 
lim g(h,(B)) J ¥(t)dt =e. 
B-o© A 


Define 


h = min [max (A, /), /,]. 
Then for h € Ha, 


A A 
(b;) F(h) < co implies h < A near co and F (h) < °° 


A 
(b.) F (h) = ce implies that Fh = o9. 
We omit the proof as it is obvious analogous of the proof of Lemma 2.14 of Jain 
and Taylor’. 


In this paper we shall denote various positive constants by the same symbol C. 
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3. Main RESULTS 
Theornm 3.1—Let 
(i) (A,) and (A,) hold, (ii) X (1) be a solution of (1.1) with EX; < ©, (iii) 
a (x, t) > O for all t and x. 
Then 
| X(t) — W(t) | = 0 (°F? (log log t)"'”) 
almost surely as f > ©9, 
Proor : We have 
> 10 decal al 0 karate: de 
t'/? (log log t)71/? t'!? (log log t)-1/? 


i a x (s), s) ds 





+ t’/? (log log) t)71/? - 


So for ¢,, = m*, A = 4/5, ma positive integer, 


eA) Wt) 1 
5d su PT os vi 
ee | ra (log log r)-1/2 |= = } 


x, 1 
7 3 su ot OX, eee ie 
ees | t'? (log log t)7?2 = am } 


j a (X(s), s) ds 


art sup Pea Mai ds Re Se v 
tmStStn 41 | aR (log log tng) 2 | > m3 
<P nea ae 
t’* (log log 1)! 2m 


fat 
tee | J a(X(s(s), sds | 
i, ee 
t |” (log log tm,1)-2!2 i 


<_ 4m* (log log tmis) ys 
tm 6 


t 
m+1 


| a (X (s), s) ds |? 


+ -4m" (log log tnes) 


tm 


(equation continued on p. 1039) 
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Cm? (log log (m + 1)4/5 


1-& 
x ( pas ) 


pessiMlogdlog-(m sul j22) 
més 


(by Lemma 1) 
< og”) 
m 


oo 
Now since Ga C (log m)/m? < oo, we have by applying Borel-Cantelli Lemma, 


P{ ap | X (t) — W(t) 


peXh(f) AACE): Ve o 
tmStStm41 11’? (log log t)~1/* ort 0) se 


and consequently 


P { lim | X (t) a) St) 


Aim “l= 731 Jog log ty | = OF =I 


Theorem 3.2—Let 


(i) (A,) and (A,) hold, (ii) X (t) bea solution of (1.1) with EX, < ee, (iii) 
a (x, t) > 0 for all t and x, (iv) h(t) > 0 increase monotonically to infinity with ¢. 
Then 
P {X(t) > P h(t)i. 0. as t > cc} = 0 or | according as 


(oe) 


h2 
I (h) = [Ae exp {— #o dt < cc or = ©. MER 
i 
Proor : Let us assume that 
h, (t) < h(t) < A, (t) for all ¢ sufficiently large asl 3,2) 


where h, (t) = (log log t)'/? and h; (t) = 2 (log log tu: 


Let us first establish the theorem under assumption (3.2). Wewill then show 
that the theorem is true for any arbitrary increasing nonnegative function h (t). 


By Theorem 3.1, for any 8 > 0 
| X (t) — W(t) | < Bt?/? (log log t)"?/ a. s. as ft > ©. 


W (t) — 8t?/? (log log t)/2 < X(t) < W(t) + BHP? (log log #)>*” 


almost surely as f > °°. Rito) 
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Let us first consider the case when 


X(t) < W(t) + Bt?” (log log t)-1/? 
<W (t) + 2 Bt*/? h- (t) (by relation (3.2)), 


Therefore 
P{X (t) > UP? A(L)i. o. as t > 00} 
se oo}, 3.4) 
qPWit) > HP (h(t) — h(t) "0 a8t> } ( 
i oe is also increasing 
Since for A (t) increasing, h (t) — Kd) : 
28 


I(h) < ce > T(h — ) < co, 


h 
So by Kolmogorov’s test criterion for Brownian motion, if / (h) < ©, then 
P{W) > fF (h@) — Fay) i: 0.08 > 09} = 0, 
Therefore, 
P{X(t)> PPP A(t)i.o. ast—o co} = 0, 
Next let us consider the case when 
W (t) — Bt??? (log log t)-2/2 < X (t). 
Since a > 0, by Theorem 3.1 the above expression can be Written as 
X(t) << W(t) + Bre (log log t)-2/2, 
So when J (h) < co, wehave _ 
P{X(t)> PR AC)i. 0. ast > co} = 0 
as shown above. 
The fact that /(h) = co = 
P{X(t)> OPA (thi.o.ast—> co} = ] 


is trivial in view of Kolmogorov’s test criterion for Brownian motion and the assump- 
tion that a (x, ft) > 0 for all ¢ and Xe 


Now let us remove the restrictj 
nonnegative function. 


Define, 


on (3.2) and consider h (t) an arbitrary increasing 


H (0) = min imax (40), b, (), be @) ...(3.5) 
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Then by Lemma 3, 
A A 
I (h) < co = 7 (h) < co and A <A near infinity. 


Again we have /; (t) < h (t) thd PS YR 
Therefore, 
Pax (tf) > PRA (iL oO. as t + co} = 0 
when I (h) < c¢. 
But h(t) < A(t) near infinity. 
Hence, 
P{X(t)>f/ h(t)i. 0. ast > cof = 0 
when J (h) < oo, 


Next again by Lemma 3, 
A 
So 


A 
P{X(t) >t)" h(t)i. 0. as as t > co} = 1. 


This implies that there exists a sequence {f,} t oo such that 


A 
bi (Be 1! * h (t,) a. 8. for every positive integer 7. 


Since I (h,) < ©°, we have 
P{X(t) > #/? hy (t)i. 0. as t > 09} = 0. 
So for {t,} t ©, 


X (ta) <& he h, (t,) a. 8. for every positive integer 7. 


Now from (3.7) and (3.8), we get 


A . 
h (tn) < hy (t,) for every positive integer n and 


hence, 
h (t,) = max [h (t,), 1, (tn)], by (3.5) 


i. €. h (t,) > h (t,) for every positive integer ”. 
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...(3.6) 


Ga) 


...(3.8) 
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Therefore by (3.7), 

dG (14 hee h (t,) a. s. for every positive integer n. 
Hence for /(h) = 0, 


P {X (t) > #2 h(t) i. o. as t > co} & 1, es) 


From (3.6) and (3.9), it is evident that we have removed the restriction (3.2). 
Hence without any loss of generality we can assume (3.2), (for the proof of this state- 
ment we have followed the technique adopted by Jain et al.® 
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ORDER LEVEL INVENTORY SYSTEM WITH POWER 
DEMAND PATTERN FOR ITEMS WITH VARIABLE 
RATE OF DETERIORATION 


T. K. Datta AND A. K. Pac 
Department of Mathematics, Jadavpur University, Calcutta 700032 


(Received 24 July 1987; after revision 27 November 1987) 


The present paper deals with a power demand pattern inventory model with 
variable rate of deterioration. Both deterministic and probabilistic demands 
have been considered. Ultimately, some particular cases regarding the de- 
mand pattern have also been discussed. 


|. INTRODUCTION 


The effect of deterioration is very important in many inventory systems. Dete- 
rioration is defined as decay or damage such that the item can not be used for its 
original purpose. Food items, drugs, pharmaceuticals and radioactive substances are 
examples of items in which sufficient deterioration can take place during the normal 
storage period of the units and consequently this loss must be taken into account when 
analyzing the system. Efforts in analyzing mathematical models of inventory in which 
a constant or variable proportion of the on-hand inventory deteriorates with time have 
been undertaken by Ghare and Schrader’, Goel and Aggarwal’, Covert and Philip’, 
Shah’, Misra® etc. to name only a few. Covert and Philip in their paper have developed 
an economic order quantity model for items with variable rate of deterioration by ass- 
uming Weibull density function for the time of deterioration and a constant demand 
rate. 


In the present paper attempts have been made to investigate an EOQ model assu- 
ming the existance of a suitable power demand pattern and a special form of Weibull 
density function. Such a special form is chosen in order to make the problem mathe- 
matically tractable. Deterministic as well as probabilistic cases of demands are consi- 
dered allowing shortages. Ultimately, some particular cases of the power demand pat- 
tern have been discussed. 


2. THE MATHEMATICAL MODEL 


The model is developed under the following assumptions : 


(i) Replenishment size is constant and the replenishment rate is infinite; 


(ii) lead time is zero; 
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(iii) T is the fixed length of each production cycle; 
(iv) C, is the inventory holding cost per unit per unit time; 
(v) C, is the shortage cost per unit per unit time; 
(vi) Cy is the cost of each deteriorated unit; 
(vii) shortages are allowed and fully backlogged; 
(viii) a variable fraction 6 (t) of the on-hand inventory deteriorates per unit time. 
In the present model, the function 6 (t) is assumed in the form 
O(t) = %t;0<a%<1,t>0 


which is a special form of two parameter Weibull function considered by Covert and 
Philip’; 


1jn ‘ 
(ix) the demand upto time ¢ is assumed to be a( z=) ’ vide, Naddor® where 


dis the demand size during the fixed cycle time T and n (0 < n < co) is the pattern 


index. (dt0-”)/")/nT1/" is the demand rate at time ¢. Such pattern in the demand rate 
is called power demand pattern. 


3. DETERMINISTIC DEMAND 


Let Q be the total amount of inventory produced or purchased at the beginning 
of each period and after fulfilling backorders let us assume we getan amount S (S > 0) 
as initial inventory. Let d be the demand during period 7. Inventory level gradually 
diminishes during time period (0, t,), t, < T due to the reasons of market demand and 
deterioration of the items and ultimately falls to zero at time f = t). Shortages occur 
during time period (t;, T) which are fully backlogged. Let /(t) be the on-hand inven- 
tory at any time ft. The differentia] equations which the on-hand inventory J (t) must 
satisfy in two different parts of the cycle time 7 are the following : 





dl t dt (-n)jn 
oe 6(t) 7(t) = — Tila» 9 Sten »-0(3.1) 
and 
dI(t) dta—nyin | 
dt = “nTime? tat ‘17: re (2 P94 | 


Solutions of the above differential equations are 


d exp ( =< es t 
tUrn) in 
nJiln 





I(t) = Sexp(— oe ys )- 
0 


x exp ( 72 ) ar, Ong fey. toea 
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and 


I(t)= Tin = ( ste = rin) igteent <7. vou( 3,4) 


Since J (t;) = 0, we find neglecting higher order terms of 60 (< 1) the following : 


fis 
S= : : bod od he 
Tin +5 Qn+hTrm "1 -45,.) 


Hence the total amount of deteriorated units 


t 
z 1 domain y dtl" Bo pe haath 
reise | nT)!" => Tin 22n+)T'™ PR ER) 


Average total cost per unit time is given by 


4 


ee Bo dt; (14-2n)/n 1 Y ie 
Cc 
C(S, ti) = ITO Ont) + T| I(t) dt — ae [ T(t) dt. 


t 
1 


Now substituting the values for / (t) given by eqns (3.3), (3.4), eliminating S us- 
ing eqn. (3.5) and then on integration we find 





a C; 60d (l+2n)/n Cid (1+n)/n 
CO) = Fan 4 NTomE 4 + Tm 
8 C Oo d pordniin = C, 8 d pOt3nln 
2 (2n aa 1) TUtn)in “1 6TUtN) jn 61 
=e C, nd plitndln Bi & C, n’ A d p(it3n)/n 
CPE DY ied Lae (2n+1) Ga+DTOmMm 1 
C, nd C, d 1/n C2 d (ltn)/n « 
tae 7 pe ft pyre 
[neglecting terms containing higher powers of 0p]. a Eg 


For minimum, the necessary condition is 


dC dC (t,) 





etatiat = 0, 
This gives 
eae 
C1 vos p+ Gi(at+ Ud, , Cro Gn+ I)d 1s 
TOtn,|n Bir n : 2n (2n +1) 1 


(equation continued on p. 1046) 
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Cin bod 3 _ C, 69 (3n + 1) d ‘oie Coal 





peas amine Per ar 6n 1 n 
Gan 
2 = T1y 

+ 24, | 


But as ¢; > O we find from the above equation the following cubic in f;. 





Le +M2+N,+P=0 ...(3.8) 
where 
Fees Cua a ] 
3n | 
| 
Cs 60 d 
a 2, os | 
N= = (C, + G), | 
‘ 
Cear | 
— ——__. | 
‘ n bie estas) 


Since (3.8) is a cubic equation in ¢, with last term P negative, it has at least one posi- 
tive root. Again since 0 < 6) < 1, we find LN — M? > 0 and so other two roots of 
(3.8) are imaginary. Let ty) be the positive root of (3.8). 


.. Optimum value of 1, is t* = to. Substituting it in (3.5), the optimum S is 


dt!" 
i 
eae Bod *(1+2n)In 
S*= Tr + Zosepre Ont Tr 4 19) 


The relation connecting S and Q is the following 


ln 
dt, 


ON de ee : (3.11) 


Hence the optimum value for OQ is 


~~ 69 d #(1+2n)/n- 
Cea S 2(Qn+ hrm ' sth 12) 


Minimum value of C is C ( t* ) 


If there be no deterioration, then 9, = 0. 


’. From (3.8) [using (3.9)] 
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er. ORE 

a OF Wie ak OF 
and corresponding expressions for Q* and S* can be obtained by substituting 9, = 0 in 
the expressions (3.12) and (3.10) respectively. 


4. PROBABILISTIC DEMAND 


In this case, it is assumed that the demand during the period (0, T) is a random 
variable x with probability density function f(x) (0 < x < co) and the demand follows 
power demand pattern with the demand rate (xt('-”/*)/nT1/". Solution of the problem 
has been investigated under the following two cases. 


Case 1: When no shortages occur 


Let us assume the inventory level of the system at any time ¢ (0 < t < T) to be 
I,x (t). Hence the differential equation which would govern the system would be 


dl.x (t) xti—n)/s 
are A tl,x (t) —— at ane”? 0 = t < fu ... (4.1) 
Solution of the equation (4.1) is the following : 
9 x oxp( — 202 : P 
2 —n)/n 9 
Ix (t) = S exp( — = Ral } esaee 7 aaa [e )/ exp ( at Jat, 
0 
A ey ee ...(4.2) 


where S (> 0) is the expected stock on hand at the beginning after meeting backorders. 
Since there is no shortage, we have 


Ix (T) > 0 
or, 
T 
S aie |" exp ( —t? ) dt >0 
0 
or, 
x S S:, 
where 
S Teh 
en ‘ 
T 


| ti-n)/a exp ( “ett ) ar 
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The expression for S, can be simplified further by neglecting higher powers of 6). We 
ultimately obtain 


O T? (4.3 
S; ~s[i- 2@n+41) ° ( ) 


The average number of items H, (x) carried in inventory per unit time is the following 


yy 
| Ix (t) dt, x < S, 


| 


Hy, (x) T 


S6. T? nx X Oo nT? 
SG ET + Catt Gntiy’® < St (4-4) 


Average number of items that deteriorates per unit time is 


] 





D, (x) = 7 Is — x — Iz (T)) 
n 6,T 
= 4ST — n+l Xs See ...(4.5) 
Average shortage per unit time is 


Case 2: When shortages occur 


If the system carries inventory during the period (0, t,) and then shortages occur 
for the remaining period (t,, 7’) of the cycle then the inventory level J,x (t) at any in- 


stant ¢ would satisfy similar differential equations as in the deterministic case and their 
solutions would be 


x exp( — $0 ; 





6 
Tnx (t) = S exp (- ome ) -- nT hyn | tQ-n)In exp ( eo t* dt, 
0 
0S ft, ...(4.7) 
and 
T (t) na vee pila 1/n 
2% Pai Tim (t, = ip hct < ya: : --.(4.8) 


Since shortages occur, we must have 


I,x(T) < 0 
or 


x > Si, where S, is given in (4.3). 
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Again, J.. (t;) = 0. This gives 


‘ 
1 


—_— ha —n)|n Go 
Ss nTiin | Seed exp ( att ar == 0, 


Expanding the integrand in ascending powers of 6, and then integrating and neglecting 
all higher order terms in 6 we get 


ep ? S 
[ + ZQn+y TOnFy 4 ]- parte (4.9) 


Taking nth power of both sides and then simplifying we find 


ACS eres iahaa(= \ar 


or 
Pe ATG ms () -..(4.10) 
where 
2(2n + 1) 2 (2n + 1) Aen 
EE et oc ee ges hake eats 
36on he NO ( x ) I. 


Solving the above cubic equation (4.10) by Cardon’s method we find the positive real 
root as follows : 


t 2(2n +1) 1 
h=u afta cae (4.11) 


where 
=4(-—G+ V(G + 4H)] 
=p [2OntD p( SV 4 {tetvr sy" 


x 2 72 
n° 6, 


fap 322 (pete) fe 


27 ng 


Expressing u and 7, ima series in ascending powers of @) and then substituting them 


in the equation(4.11) we find 


Payee uel n0o »(=)” 12 
n= (2) uf a tr) x AGE 


neglecting all higher order terms in @). From the above equation the expressions for 


rd ig Ya , ‘eed " etc. which are needed to calculate the total average expected 


cost can be obtained as follows : 
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3n 
In S Pe alc oes asamin( >| ass( 4,13) 
6 OS Sek Ona x ( 
n S Bn+1 
Gi+n)/n S = n)/n 21g 1) 80_ cisamyn ( = ) 
bn ( ) en AB Cee ae x 
...(4.14) 
n)/n Soe n)jn 21 3n + l (1+5n) Jn (=) 
tah =(=) ee TEN on tt ee 
...(4.15) 


The following expression for S in terms of t; can be obtained from eqn. (4.9) 


I/n (i+en)/n 
t Bo t\ 





1 ; 
ag [ Tin + 2 ere sis ...(4.16) 
The average number of items H, (x) carried in inventory per unit time is the following: 


t 
1 


H, (x) = val Ix (t) dt, x > S, 
0 
= x I G+n)/n 06 (i+3n)/n 
—- TGtnrin E + 1 ty “he 3 (3n+1) ty a x = S}. ...(4.17) 


Average amount of inventory that deteriorates per unit time is 


t 


1 
l xtG-n)in 

D, (x) = +(s — ary dt |x>s, 

0 


s xin 
~ [ _ Patan: ...(4.18) 


Average shortages per unit time is 





y 
1 
Gi ae | le (aS 
y 
® n 1 
(tn)/n __ I/n Gtndln 
Tatn)jn E say, TUT} Tt, + Peep tit. ik 
x > Sy 


...(4.19) 
’. Expected total cost of the system per unit time becomes [using (4.3) for S,] 
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pale ron ] 





Ki, S) = D, (x) f (x) dx + Cs 
x 7 D; (x) f (x) dx + C, 
s[1- eo Sn 
2 (2n + 1) 


eee as: ae T* 
Po at ot + 1) Tea | 
A, (x) f (x) dx 


+ C H, (x) f (x) dx 


Sor 
se om 1) 


ts Pent ] 


+C G, (x) f (x) dx 


i G, (x) f (x) dx. 
Ay aie Fa ihe 
2(2n + 1) 
Now substituting the values of D, (x), D2 (x), H, (x), Hz (x), G: (x), G2 (x) from (4.5), 


(4.18), (4.4), (4.17), (4.6), (4.19) respectively and finally eliminating t, using (4.12), 
(4.13) (4.14) (4.15) we get the following : 














6 T? 
s[1- 2 (2n + 1) 
K (S), = CG, [ 3 so 7 — ix |p (x) ax 
0 
co 
21 A on 
+ Cs | hs mites | 
6) T? 
s[1— 2 (2n + 1) 
x f (x) dx 
0 T? 
Sloan 
SO, T? nx 
he | [s- 6 atl 
69 n> T? ] d 
+ G+ DGrth* 7” 


(equation continued on p. 1052) 
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ii 1- sett 
+ C, | = 
0 T? 
s[1- 2 (2n + 1) 





(157n + 47) 9% Pa Sok | d 
" 98Gn + 1)Gae 1) one eee 


co 
nx 218 f he 
te fesse 16 (2n +1) 
60 T? 
2 [1 2 (2n+1) ] 
se | snt 216) 7? Sett_ ] : 
ait Tet. 16a al) Gee ee 
...(4.20) 


If the probability density function / (x) and pattern index n are prescribed, then right- 
hand side of eqn. (4.20) can be evaluated. The necessory condition for the minimum 
expected cost K (S) is the relation 


dK (S) _ 0 
DSF ae 
Equating ds ‘0 zero, the optimum value of S = S* (> 0) can be derived. For 
ive d? 
this value of S = S*, the sufficient condition for minimum ES st >0 would 


also be satisfied. 


5. Discussion 


In the present problem, a power demand pattern has been assumed with demand 
rate (dt\-™)!")inT"/", where T, d, t are prescribed cycle time, entire demand during 
(0, T) period, ¢ is time (0 <t < T) respectively and n is pattern index. Substituting dif- 
ferent values to the pattern index 7 in the equations for the total cost and total expect- 
ed cost per unit time given by the equations (3.7) and (4.20) respectively, we can deter- 
mine the corresponding cost equations. Then differentiating and proceeding in the 


usual manner the optimum values S*, t* Q* etc., can be evaluated. Substituting n = 


and n = co in the power demand pattern formula it can be seen that these two cor- 
respond to the two extreme cases, i. e., when the entire demand occurs at the end of 
the period and when the demand is instantaneous in nature. Forn = 1 the demand 
pattern is uniform and if there is no deterioration then it corresponds to the case dis- 


cussed in Wilson’s model. In this case the total cost per unit time (deterministic case) 
given by (3.7) reduces to [by substituting n = 1, 0 = 0] 


C; d Cid d 


C(t) = co a sis + (C, + C,) WEE Crd 


fi 





ti — th. 
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dC 


Fo i —— 
r optimum C, dt, 


= 0, 


Equating the derivative to zero and Simplifying we find the optimum /, as 





CT 
ety ee 
, : C,+C, 
f d?C Pe 
ince > 0, C would be minimized for 1, = ¢* . 
Ct ed a 2 


For n = 1, 60 = 0 we find from equation (3.10) 


dt* 
s* = van G Ce ee i 
is Te ete 8 iG a 


Similarly, the expected total cost per unit-time (probabilistic case) given by (4.20) re- 
duces to 


s 


K(S) =C, | [s a > |F@ ae eS 5 P 
0 S 


sive Wee ss _ J 100 ax. 


Ss 


Solvin ae (Sy = 0 we can determine the value of S*. 


When n > 1, a larger portion of the demand occurs towards the beginning of the period 
and when 0 < n < 1, a larger portion of the demand occurs at the end of the period. 
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We consider the differential equation 
Lpl=- hee +g (X)y() =2Y (XI), 0S x < 
° 


where 4 is a complex parameter. It has been shown that the convergence of 
the expansion of any function f(.), square integrable on [0, oo), in terms of 
the eigenfunctions of a boundary value problem associated with the differen- 
tial equation is independent of the function g (.), provided g (.)is square 
integrable on [0, 00). 


§ 1. We consider the differential equation 


d’y (x) 


L [y] os dx? 





+ q(x) y(x) =Ay(x),0< x < © Late} 


where A = u + ivisa complex parameter. One of the usual problems with such a 
differential equation is to consider the convergence of the expansion of an arbitary 
function f (.), square-integrable on [0, oo) , [we write f © L* [0, co)], in terms of the 
eigenfunctions of boundary value problem associated with the differential equation 
(1.1). The idea here is to prove the following : 


Theorem—If q (.) € L? [0, ©), and the eigenfunction expansion of an arbitrary 


f € L*[0, cc) associated with a boundary value problem’ (I]) consisting of the differ- 
ential equation 


; d® y (x) 
by DV] me Ay eee mm tbe) 


and the boundary condition 


y (0) cos « + y’ (0) sinn = 0 O<a<n) cai 123) 


is convergent, then so is the eigenfunction ex 
value problem (I) consisting of (1.1) and the 
boundary value problem (II). 


pansion of f associated with the boundary 
same boundary condition (1.3) as in the 
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To establish this result, we obtain a relation between the ®-functions (for the 
definition of the ®-function we refer to section (2.6) of Titchmarsh’ for these two 
boundary value problems and show that the contribution corresponding to the 
difference of these two ®-functions towards the expansion of / is nil. 


§ 2. Let 6 = ¢ (x, A) and 6 = 6 (x, A) be the two linearly independent solutions 


of the differential equation (1.1) satisfying the following initial conditions : 


¢ (0, A) = sina, @ (0, A) = cosa 


$9) (0, A) = — cos a, 90) (0, 2) = sin « —nan7<a<cn, WAAL 


It has been shown by Titchmash’ (§ 2.1) that there exists a function m (.), analytic in 
the two half planes im A > 0 and im A < 0 such that, for im A + 0 
b (x, A) = 8 (x, A) + m (A) $ (x, A) E L?[0, ©9). se ¢y)) 


With the help of these two solutions ¢ (x, A) and ¢ (x, A) for any f € L?[0, co) let us 
define the function ® (x, A; f) = ® (x) as 


® (x) = d(x, a) f $n, DF) dy + 8% 2) F VOLO) a. 


This is the so called ®-function of Titchmarsh associated with the BVP (]). It can be 
easily proved that ® (x) satisfies the following differential equation 


L [y] = Ay (x) —f (*). pr(2-3) 
Titchmarsh’ (§ 3.1) has proved that the expansion of the arbitrary function 


f € L? (0, ce) associated with the boundary value problem (I) is obtained from 


R,is 
pl | rdhix, Ay f id A ...(2.4) 
‘ —R+1& 


by making R > co andé > 0. 


Let do (x, A), Oo (x, A) be the solutions of the differential equation (1.2) satisfying the 


following : 
. Q) 
= 0,A) = — cosa 
fo (0, A) sina, ¢, (0, A) PAN PSE, 
6 (0, A) = cose, me (0, A) = sin « 


It follows easily that 


| bay: A 
8) (x, A) = cos a cos (x V/A) + ae sin « sin (x, V/A) 


$y (x, A) = sin a cos (x VA) — x cos « sin (x V/A). 42.5) 
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So, the square integrable solution ty) (x, A) of this differential equation is given by 
(Titchmarsh’, § 4.1) 


cix/a er 2.6) 
ty (x, A) = iE EEA ( 


As q € L* [0, 0°), we know that the differential equation (1.1) is separated i.e. f and 
L[f] € L? [0, cc) should imply that gf € L? [0, oo). 


Since for each f € L? [0, co) we know that ® (x, A; f) € L* [0, ce) it is follows that 
® (x, A; f) and L [® (x, A; f)] = A (x, A; f) — f E L*[0, co) and by seperability of 
the differential equation (1.1) we have q © (x, A; f) € L? [0, c0). Now 


® (x, 2, g © (x)) = Yo (% AF $0 (A) g (9) © (9) dy 


+ $0 (x, a) F do (v, 4) g (A) © (y) dy. Can 


Since ® is a solution of the differential equation (2.3), we get 
L [® (x)] = A ® (x) — f(x) 
or 
Ly [® (x)] + 9 (x) © (x) = AO (x) — f(x) 
or 
Lo [® (x)] = A ® (x) — g (x) 
where g(x) = f(x) + q(x) ® (x) € L?[0, cc). This implies that © (x) is the 


®-function of the boundary value problem associated with the differentia] equation 


(1.2) corresponding to g (Note that both the boundary value problems have the same 
boundary conditions). So, — 


® (x, Asp) a Dy (x, A; g) ap ®, (x Af + q ®) 
= x; Az) + , (x, A; q ®), wa0(2.0 


In view of (2.4) and (2.8), it now follows clearly that in order to establish the theorem 
we need only to show that the contribution to 


wards eigenfunction expansion arising 
from Dp (x, A; gq ®) is nil, 


The following two lemmas will be useful for further discussion. 


Lemma 1—For any f € L? [0, co), A = y + iv, 


J low@anirdr< 4 | 1s) | eax. 
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For proof we refer to § 2.8 of Titchmarsh’. 
Lemma 2—If WA = s = 6 + it,A = u + iv then 
®, (x, A; g D) = o( sat): 
lvl |s| 


PROOF: 


“« 


ixs Ps l s 
c | {sin “ cos (ys) — 5 COs « sin os} 
0 


® MCE LO) ese = aie a a 
alte 2) cos « + is sin « 


q (y) ® (y) dy 





co 
; | ‘ 
4: { sin « cos (xs) — 5 00s % sin (xy) ety, 
cos « + is sin « - 
; ] : 
Now sin « cos ys — 5 008 & sin ys 
sina, , Cos & : 
(6 |e 7?) — 2 (e'?*, — .e- 1! 














\| 


Sin « _ COS & iys (3 ae) —tys 
(3 e+ ae ora ke 


























P _jyy | Sim @ — cos« ty [sina , cose 
=O (ew 2 S|) +o (e 2 |) 
= O (e!”) 
So, 
; ; COs K-. } 
e**! + sin « cos ys — sin ys 
= O (e-'*) O (e’’) = O(1) sineeO Sy Kx. 
Similarly 
sin « cos (xs) — ao gin (xs) = O (e!*) 
and so 
{sin % cos (xs) — = * sin (xs)} e1t = QO (e* e~*) for all y > x 
= O(1). 
Hence 
1 oo 
0 (24590) = O( et ey | 120) 014) ) 


0 


(equation continued on p. 1058) 


1058 JYOTI DAS AND ANINDITA CHATTERJEE 


co 


of L{fiaorirah tl 1, as) 12 av} | 


= O ( TF ), using Lemma 1], for all A. ...(2.9) 


§ 3. Now consider the contour I (R) formed by the segments of lines (— R + i, 
—R + id),(—R + id, R + i5), (R + i5, R + i) joined by semicircles of radius R and 
centres + i. Using (2.8) we get, 


D(y, Sf) dA = § D(x, rArIf)dA — f O(x,3q0)d. — ...B.1) 
I(R) T(R) r(R) 


Using Lemma 2 we get 


0, 2g) OO) — o ( | SENGEH dn ). 3.2) 


T(R) (R) 


On the part of the upper semicircle in the Ist quadrant, we have 


A= i+ Refs (0<¢< =) 


Hence, | ®, (x, A; q (x) © (x, A; f)) dA integrated round this semicircle in the first 


quadrant, gives, 





_ ™12 
Rd @ Z Rdg ) _ 
‘ {| RP (Ll + Rana we 1 TE Rane} = 0(l) as R > co, 


0 
Hence finally it is a question of proving that 
7 R+l 
lim iB ®y (x, A, g (x) D(x, A; f)) dA = 0. me RZ) 


830 R+i 
R-.©o 


Using Lemma 2 we get 


Ril R+i 

J 0 (x As 4 (x) © (x)) aa = | of 43 )a 

ah lVvaAl |v], 
R+is 


1 


= 


O dv ; : 
re [R? +. y2}t 14 | y | ), since A =u + vi 
& 


(equation contiuned on p. 1059) 
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1 
= dy 
JO paar? 
§ 


(R"41)'/4 
at) ( | eet aii 2S 
(t? —R) (12+ R) ), 


(R°4.82)1/8 


substituting R? + vy? = ¢4 





ay : 
= (t? — R) + (t? + R) 
° ( — R) (+R) dt ) 
(R'48°)1/4 
= atts Sy Ate ces IE oe (RE +. Bt A 
azaks VR iin meno 
+0 (+ pga Etat ah SBC a a v)) 
VR 18 CREE IE VR RE BER} 
1/4 52 \1/4 
ae Ca) 
me Oe tay te th 
VR 1 \1/4 82 \1/4 
1 + (1 + wr) (3 ae =) 
1/6 52 \1/4 
ai a —1 ee Se a a | 
u (1 as ( R ) 
+0 JR log 


I \'A ss] gst toate: 
(1+ a aa (1+ 5) “4.9 


= 0(1) as R > oo; 3 + 0 by choosing 3 suitably as dependent on R. 
Hence (3.3) follows. 
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CLASSICAL LIE SUPERALGEBRAS 4 (m,n), B(m,n) AND D (m, n) 
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Modified Satake diagrams are obtained for the Lie superalgebras A (m, n), 
B(m, n) and D(m,n) in order to get the involutive automorphisms of these 
Lie superalgebras. Iwasawa and Langlands decompositions are done on the 
same lines as for the ordinary Lie algebras. These decompositions facilitate 
the determination of the parabolic subalgebras. The parabolic subalgebras 
can be used to determine the representations using Schmidt constructions. 


1. INTRODUCTION 


Graded Lie algebras or Lie superalgebras have generated a great interest in par- 
ticle physics in the contex of supersymmetries. In this paper, we make a preliminary 
attempt to obtain the induced representations of some Lie superalgebras by the Schmidt 
construction’ method. In order to do this, the minimal parabolic subalgebras and 
other parabolic subalgebras have been found out through which one can get the Lie 
Supergroups corresponding to some Lie superalgebras by induced method. Very little 
is known on the representations of graded Lie groups except Kostant’s novel method?. 
Attempts have been made to obtain typical and atypical representations of some special 
Lie superalgebras. The analysis of irreducible representations seem to be inadequate®. 
Here, we attempt to present.a systematic study of the representations of superalgebras 
through the prescription subscribed for ordinary simple Lie algebras. 


We devise here the super Satake diagrams for the Lie superalgebras on the same 
lines as for ordinary Lie algebras*’®, to find out the automorphisms of the Lie superal- 
gebras A (m,n), B(m, n) and D (m,n). The plan of the paper is as follows. In Section 
2 following Kac*”’, we give a brief and quick resume of the classification of classical Lie 
superalgebras along with the root systems and Dynkin diagrams of Lie superalgebras. 
In section 3, we give the real forms of classical Lie superalgebras. In section 4, we 
recapitulate the salient features of Iwasawa and Langlands decompositions for Lie su- 
peralgebras in the light of earlier analysis for classical Lie algebras®-?°. In section 5, we 
discuss the modified Satake diagrams for classical Lie superalgebras A (m,n), B (m, n) 
and D (m, n) of physical interest. In section 6, we display the Iwasawa and Lang- 
lands decompositions and obtain the parabolic and minimal parabolic subalgebras. 
From the knowledge of parabolic subalgebras, one can resort to Schmidt construction 
and derive the corresponding induced representations for superalgebras. 
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2. CLASSIFICATION OF SIMPLE CLASSICAL Lig SUPERALGEBRAS 


Here, in this section, we present a quick survey of finite-dimensional classical 


simple Lie superalgebras over C (complex field), their root systems and Dynkin diagr- 
ams analysed by Kac. 


Let Vo = V @ Vz be a z,-graded Space, dim Vo = m, dim Vz = n. Let L(V) or 

1 (m, n) be the vector space of all the (m +n) x (m + n) matrices, written in block 
B ; ; ; 

form X = ( CD ) where 4 Is an arbitrary m x m matrix, Ban arbitrary m x n 

matrix, C an arbitrary n ~ m matrix and D an arbitrary n <n matrix. The Lie algebra 


8 Of / (V) consists of the diagonal block matrices ( . ) The odd subspace gr 


consists of off-diagonal block matrices A O ) The bracket < x, x > is commut- 


ator of two elements of / (V) if x or x’ is an element of gz and its an anticommutator if 
x and x are elements of gz. 


(a) The Lie Superalgebras A (m, n) 


The special linear graded Lie algebra s/(m, n) consists of block matrices ( A 


such that TrA = TrD, which is a graded ideal of / (m,n) of codimension one. Its Lie 
algebra is s] (m) x sl (n) X gl(\), (g/ (1) is the trivial one-dimensional Lie algebra). 
If n + m then s/ (m, n) is simple. 

We get 


A (m,n) = sl(m+1,n+ 1) form#~n,m,n> 0. eet) 


The roots of A (m, n) are expressed in terms of linear functions ¢,, ..., €ms1, 91 = Emaos 
--+» 8mt1 = €mtn+2- Even roots are given by 4, = {e, — €,; 8, — 8,},i 4 j and the odd 
roots are A, = + (¢ — «,). The simplest system of roots is 


{e, — &, €g — €y, .00y EmtTy s,— 3, a Se, Boag Sn41}- mA ty) 
(b) Orthosymplectic Lie superalgebras 


Suppose that m = 2r is an even positive integer. Let the subalgebra of / (m, 2r) 


A B : : 
consist of all block matrices ( CD ) which satisfy 


taG4GA-0, tp+Dao, C = tao. et (aee) 


This subalgebra is simple and its Lie algebra is Sp (2r) x O(m). It is denoted by 
osp (m, 2r). 


(i) Lie superalgebra B(m, n)—By using Cartan’s notation 
B(m, n) = Osp(2m + 1, 2n),m>0,n > 0. ...(2.4) 


The roots of B (m, n) are expressed in terms of linear functions €, .. , €n, 9) = €gm415 
.--y On = €om+n. Lhen even roots Ay are given as 
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A,g={zateptBWitagtkat Spi FT; 
and the odd roots are 

A, = {+ 83 + + €)}. 
The simplest system of simple roots is given by 


{81 Bs 8e, soe9 8, — €1, €) — €gy .-. €m—-1 — €ms Em} ifm > 0, 
and 

{3, a Se, eres Sy-1 a ai Sn} A Ae) 
ifm = 0. 


(ii) Lie superalgebras D (m,n), C (n)—The Lie algebra of D (m,n) is of the 
type Dm ® C, and it can be written as 


D (m,n) = Osp (2m, 2n), m > 2,n > 0. Rae 


The roots of D (m, n) are expressed in terms of linear functions €,, ..., €m, 51 = €2m415 
weep On = €oman. Even roots are 


A={taete; $2; 4+545},14/ 
and the odd roots are A, = {+ « + €,}. The simplest systems of simple roots are 
{85 = Oey sis Op egy Cas ee eee Em-1 + Em} 
(0; i Sassen ae eg ae en ee, OE Rieke 
The Lie superalgebra C (n) is defined as follows. 


The matrices of the Lie algebra of C (n) are of the form 


m 0 





— 





a b 


eerie, 
ee RE ee 





c —a' 


where a, b and c are (n — 1) x (n — 1) matrices, b and c being symmetric and « € k, 
C (n) = Osp (2, 2n — 2),n > 2. ...(2.8) 


The roots of C (n) are expressed in terms of linear functions €,, 5 


1 = €3, «.-5 5y- = €n+1- 
Even roots are A, = 


{+ 28); + 3) + 8)}; and the odd roots are A, = {+ «, + 3}. 
Following are the systems of simple roots 


+ {e, aa 81, 8, aoa on tery Sno» _ Gey. pda & 
a {3, a be, wavy OF — Si, Sf S141, sory Ons 7 Ost, 28,-1}5 


== {8, i So, tery 8,-2 = Bors eer me Si 8-1 + «,}. .-.(2.9) 
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(c) Lie superalgebras QO (n) , n > 2— Define a subalgebra O, of sl (n+ 1, 
. B = 
n-+1) by matrices of the form ( BA where 7r B= 0. The quotient algebra Q,/z 


Vv 
= Q (n) is simple where z is the centre of Q,. 


(d) Lie superalgebra P (n),n > 2 —This is also a subalgebra of s/ (n + 1, n+1) 


containing matrices of the form ) where 7rA = 0, Bis asymmetric matrix 


C —A! 
and C is a skew-symmetric matrix. Its Lie algebra is s/ (n). 


(€) The Lie superalgebras F (4), G (3) and D (2, 1; «)—There is one and only 
one 40-dimensional classical Lie superalgebra F (4) for which F (4) is a Lie algebra of 
type B; ® A, and its representation on F (4) is spiny @ s/,. The roots are expressed in 
terms of linear functions ¢€,, €2,€; corresponding to B; and 5, corresponding to A). 
Then even roots are Ay = {+ ¢ + €;; + «3 + 5}, 


i ~ j; and odd roots are 4, =4 (+ €, +, + €, + 3), 
There are four systems of simple roots : 
Silent pone ag at) Gi 4g; <,); 
{—5, (44 +¢« +6, + 3), — &,€, — €}; 
tehtiel  Giee *), + (— ey + + oo — 9), 4(— 4 — Gg — & +9), 
€; —.€5]; 
{h(a +e, + €; +), $ (a: — € — es — 9), €2 — 1, €3 — &}. —...(2.10) 


There is one and only one 31-dimensional classical Lie superalgebra G (3) for 
which G (3), is a Lie algebra of type G, ® A, and its representation on G (3)z is Ge ® 
SI,. The roots are expressed in terms of linear functions €,, €2, €s, corresponding to Gs, 
«, + €, + €; = 0, and 4, corresponding to A). 


The even roots 

Ay = {a — €;; + € + 23} 
and the odd roots are 

A, = {+ « + 8; + 3}. 
There is a unique system of simple roots 

{5 + €1, €2, €, — €2}. 


There is one parameter family of 17-dimensional Lie superalgebras D (2, 1; «), 
a € k*/{0, — 1}, consisting of all simple Lie superalgebras for which D (2, 1; «)g isa 
Lie algebra of type A, ® 4, ® A, and its representation on D (2, 1; «)z is sl, x sl, 
@ sl,. The roots of D (2, 1; a) are expressed in terms of linear functions 41, ¢, and €3. 


(2.11) 
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Even roots are A, = {+ 2 «} and odd roots are {+ «1 + €2 = es}. There are four 
systems of simple roots : 

{e, + 2 + €3, — 2, — 2€,},1 A J; i, j = 1, 2, 3; 

{e, + es + €s, €; — €y — Cg, — 41 — My F As}. ) ee A 


After classifying all the classical simple Lie superalgebras, we tabulate the Dynkin 
diagrams of these superalgebras in Table I. The O, @ and @are called white, grey 
and black respectively. The white circles imply even roots whereas the grey and 


TABLE I 
Kac-Dynkin diagrams of Lie superalgebras 


A (m,n) O——0-----8——0-----O 0 
B(m,n) O——0-----@—0-----0=3=0 
B(o,n) 0 Nae 9 ae oi foes 
C(n) ee Wee ee ka ee 
D(m,n) eee ee a gota 

D(2,1,@) 

F{4) o—_O=S0—_O | 

G(3) oie 


black ci 
i Saree ris ve roots, The roots are expressed in terms of linear functions e;, 
ale ie sabi a» 009 Oy which form a unit basis of h*, the dual space of the Cartan sub- 
i<mi<k ae Product (st, ¢)) = Bi, (84 8) = — Buy (e;, 5.) = 0 where 1 < i 
S ’ Cy < ; ; = a ? ’ ake ~~ * 
= Sms If he G ntl 2) ree rank of the superalgebra) are the 
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generators of the Cartan subalgebra and ar ( a ) are the generators corresponding 
to the ith positive (negative) simple root, then 


E Pax | = bi, hy, [hu, hy] = 0, ai, at ] enya (2.13) 


where a; are the elements of the Cartan matrix. The remaining generators may be 
defined from the simple roots by (anti-) commutation. 


3. REALS Forms oF FinttE-DIMENSIONAL SIMPLE CLASSICAL LIE SUPERALGEBRAS?’. 
Let g = gy ® gi bea complex classical Lie superalgebra. A real Lie super- 

algebra ge of g is a real form of g it g is the complexification of ge. Such a real form 

ge determines a mappings: g—g. This mapping o has the following properties : 


1. is semilinear i.e. [¢X, oY] = o[X, Y] for ¥, Y € g. 
2. ois an involution i.e. co? = /g. 

The involutive automorphism c can be written as o = 6, + 6, where oo, 5) are 
its restrictions to gy and gz respectively. This can be shown that Zo = Loe DB Bre 
where gge = {x + ox | x E gg} and gre = {x +ox| x E gr. Wealso see that 
ifs and o’ are two involutive automorphisms of g then the real forms ge and go’ are 


isomorphic iff there exists an automorphism ¢ of g such that o’ = ¢o¢-. Also every 
inner automorphism of gz extends to an inner automorphism of g. 


The real forms of the classical Lie superalgebras can be uniquely determined upto 
an isomorphism, by the real form gye of the Lie subalgebra g>. The real forms are 
listed in Table II from where one can single out the compact real forms of classical Lie 
superalgebras. 


4. IWASAWA AND LANGLANDS DECOMPOSITIONS OF SIMPLE CLASSICAL LIE 
SUPERALGEBRAS 


Let g = gy ® gr bea complex classical Lie superalgebra. Then there is a com- 
pact real form ge of g such that 


g-k@ip, g=k @p, we Ce) 
where k = (1 + sc) ge and p = i(1 — ¢) ge, and ¢ is the involutive automorphism. 


Let h be the Cartan subalgebra of gz and A the set of positive even and odd 
roots, (Az U Az) of g with respect toh. Following commutation (anti) relations are 
satisfied by the elements of g, 


[e., kh] = a(h)ex., hE handa € A. 
[ex, en3] + 0, if«, Banda +BEA 


= 0, otherwise. 
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TABLE II 


Real forms of Finite-dimentional classical Lie superalgbra 


nn UT ETEETSESSESSEESSISEIESSSSSS am 


Lie Superalgebra g Real Form 8o0¢ 
A (m, n) 1. Su* (m) @® Su*(n)@ R 
2. Su(p,m—p) ® Su (r,n—r) @ iR 
3. sl (m,R) @® Sl(n,R) OR 
4. sl (n, C), ifm=n 


B(m,n), D (m, n) Sp (m, R) ® SO (s,q) 
Sp (r,s) ® SO* (2p) 
Sp (n, R) ® SO (2) 

2. Sp (r,s) ® SO (2) 


Sp (m, R) 


WN US 


C (n) 


— 
. 


B (m, 0) 


— 
. 


P (n) Su (n), if nis even 


SI (n, R) 


nN 


Q (n) 1. Su(p, n—p) 
Su* (n), if nis even 


3. Si(n. R) 


v 


D (2, 1; «) 1. sl(2,R) @ si(2, R) @ sl (2, R) 
2. su(2) ® su (2) @ sl (2, R) when « is real 
3. sl(2, R) ® sl(2, R) whene +g = —1 


G (3) sl(2, R) ® Gayo 


sl (2, R) ® Ga, 
S1(2, R) ® SO (7) 
SI (2, R) ® SO (3,4) 
SU (2) ® SO (2, 5) 


Oy aes 


F (4) 


YN 


[ex, ea) = (ea, en) he. 
ahs ---(4,2) 


(hag b)imite (hh eon 
We define the Killing forms as 

B (a,€s2) ey 
The compact real form ge may be taken to consist of thy, for « = 


of g together with (e. + i = a, ..., a, J = rank 
e- = 
ah 0 «) and / (¢, — e_g) for every root « of g defined with re- 
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The basis element ih., for « = « pJ = 1, ...,1 (1 = rank of the Lie superalgebra) 
correspond to eigen value + 1 while (€x + @-«) and i (e, — e_,) both correspond to 
eigen value + | (exp {x (A)} = + |). kK has a basic consisting of ih, together with 
(ex + e-«) and i (e. — e_.) for all x such that exp {« (h)} = + 1, while the basis of p 
consists of i (e. + e_x) and (e. — e_x) for all « such that exp {a (h)} = — 1, ais the 
maximal abelian subalgebra of p with dimension m,, and its basis may be taken to have 
elements of the form i (e. + e-.). Now choose m, the centralizer of a in k. Its basis 
elements are of the form (e, + e-,). The complexification of (a @ m) gives a Cartan 
subalgebra of g. This second Cartan subalgebra of g, denoted by h’ determines the 
nilpotent subalgebra n as follows. 


There exists an inner automorphism, V that maps h’ into h such that h, = Vhi 


whereV = | Vo, hy € bh andh, € hand the product I Vq is over all «’s con- 


a a 


Stituting a and m. Let A* denote the set of positive roots defined with respect to h. 
Then 
h, = 2 b; («)h,and« € A* if b, (x) > 0. ...(4.3) 


and if / is the least index, then b, («) #0. 
We divide the positive roots A+ in two classes as follows : 
Ps = {aa © At, a(h) # a (VoV-1 hy} 
p~- = {a:a€ A*,a(h) =a (VoV"h) Y¥ hE hh. ...(4.4) 


Let the subalgebra n be spanned by elements V-' e, fora € py andn =n/Q g. 
n is the nilpotent subalgebra of g. Thus the Iwasawa decomposition of g is given as 


g=k@a@n. wask4e5) 


Now we discuss, the Langlands decomposition to obtain the parabolic subalge- 
bras of the Lie superalgebras. 
A minimal parabolic subalgebra is defined to be any subalgebra that is conjugate 
to 
Pa—m@aen. ...(4.6) 


Any subalgebra of g containing a minimal parabolic subalgebra is a general parabolic ~ 
subalgebra. There exists 2!™,! conjugacy classes of parabolic subalgebras of g and in 
each such class there is a standard parabolic subalgebra which can be obtained as 
follows. 

Let = be the set of roots for a and let / be the set of positive rootsin 2. Let @ 
denote the subset of . Let <@> denote the set of roots in = which arises as linear 
combinations of roots in 6. We define 
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<9> =2,1 <O>, <9>-=2_N <> .»(4.7) 
where 5., - denote the positive and negative roots in &. 
Let n, (9), n- (0), 1 (8) denote the subspaces of a corresponding to <@>., <6>- 
and 3, — <@>,. Now let us define 
ae = {a € a, \(a) = 0, for all A € 98} ...(4.8) 


and let a (6) be the orthogonal complement of a» in a with respect to the Cartan Kill- 
ing form. Then 


Po = Moe @ ao @ No oA CD 
is a parabolic subalgebra of g, where 

me =m @ m4 (6) @n. (0) @ a (8). ...(4.10) 
A rea! Cartan subalgebra h is said to be o-invariant if 

h = (h () k) @ (bh p). .(4.11) 


A parabolic subalgebra pe is said to be cuspidal if there exists a o-invariant real Cartan 
subalgebra h such that 
Qo =h/f) p. ...(4.12) 


This shows that the minimal parabolic subalgebra is cuspidal. 


5. SATAKE DIAGRAMS! AND INNER AUTOMORPHISMS!2 OF CLASSICAL LIE 
SUPERALGEBRAS 


Let ge be real simple Lie superalgebra and g, which is a complex simple Lie su- 
peralgebra, be the complexification of ge. Let h be the Cartan subalgebra and h* be 
the vector space dual of h. A is the set of roots of g with respect to h. Let C be the 
conjugation of g defined by ge so that 


C (x1 + ix.) = x, — ix, for x,, x. © ge : ese(5.1) 
C acts on the root space A as follows: 


for each root « € A, define c («) by 


o(a)(h)= a(C(A)),hEh (S21 
then 
C (g*) = ge (*), 


The mapping « > o («) extends by linearity to an involutory isometry of the Euclidean 


ee under which A — A, is stable, and Ay is the set of roots « € A. such 
that” : 
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a = (—1)/*Ilo (a) et RD, 
| « | is equal to zero for an even root and is + 1 for an odd root. 
We have 

“a — (—1)!*! g(a) = &foralla € A. et:4) 


Let B be the basis of A such that By = BO Aj is a basis of Ay, and such A+ — AS is 


o-stable, where At and Aa are the set of positive roots determined by Band B, respec- 


tively. The involution o determines an involutory permutation of B — By as follows : if 
« € B— B, there exists a unique 8 € B — By such that 


o (a) = B, as.) 
and the mapping « + 2 isa permutation of order 2. We have 
P (a) = P (8) = } (a + 8), where p = 4(1 +), .--(5.6) 


Each real simple Lie superalgebra go determines a normal pair (A, s) which determines 
< upto isomorphism. A pair (A, ) is said to be normal if« €E A>a+oa(a) A. 
The root space A may be represented by its Dynkin diagram, the vertices of which re- 
present the elements of the basis B of A. The action of o may be indicated as follows: 
The vertices of the diagram which represents the elements of B, are coloured black and 
the remainder are coloured white. Two white vertices aes elements «, 8B € B 


— B, for which P (x) = p (3) are joined by an arrow 0 0. The resulting diagram 
is called the Satake diagrams of Be and determines gz upto ismorphic Satake diagrams 
for the Lie superalgebras A (m,n) m ¥~ n, B (m,n) and D (m, n) are given in Table III. 
The involution o of A is uniquely determined by the Satake diagram. 


TABLE III 
Satake diagrams for A (m, n), B (m ,n) and D (m, n). 








Lie Superalgebra Satake Diagram 
A(m,n), m#n a fons ae 
B(m,n) Mie oe koe 

od Xs n?7"n ntm-1 Nem 
“Sn 
o—------ @——®------- ae 
pine A, Mp a2, 
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Let 
BEB ee a, cy) 
and 
By = {B:, ..-, Bs}- eAgal? 
Then one can show that* 
—o (a) = any + 2(—1)!F lau Bi, B E Ao -te) 


where 7 is the involutive permutation of 1, 2, ..., r and my are non-negative integers. 
6. EXAMPLES 
In this section, we illustrate our analysis for superalgebras of physical interest. 


(a) (i) A (0,1)—The positive roots are o1, % and %; + %. exp {a (A)} = + 1 for 
a, and exp {a (h)} = — 1 fore, and a, + a. Thus k hasa basis ih, for « = %, a» 
and (e. + eu), i(@x — e-«) for « = a1, p has a basis i (e. + e-«) and (¢. — e-«) for 
a= ho, %, + &. a consists of element 


ht = i (€a + e_«) for a= ho ..-(6.1) 


and m has a basis 


— th, = @a + @-a, % = %. = .( Ore) 
Thus Ra = {o.} and Rm = {a,}. Therefore V = V. with « = «a, and it gives 
h, = — hy,, h, = + he,. Using eqn. (4.3) and eqn. (4.4) we get 

A* = {a;, %, + a, — a} 3..(6.3) 
and 

P, = {a;, — o, @ + a}. ...(6.4) 


The basis elements of n are given by n = n () g where the elements of n are V-’ e, 
for « € p,. We see that for this superalgebra we have only two parabolic subalgebras. 


(ii) A (1, 2)—The positive roots of A (1, 2) are a,, a, a3, a4, a1 + as, % + a, 
Ost Og, &, + Oy + ag, % + a3 + a, and a, + a + a3 + a. exp {« (A4)} = + 1 for 
X= Hy Gay Oy H Oy, & + % + a3 whereas exp {a (h)} = — 1 for a = as, a4, ata, 
3 + Hy, Xe + as + ae and a, + a» + a, + ay, Therefore, the basis of k consists of 
ihe for « = a1, a, a, a4, (Cn + €-a), i (€C1—€_a) for « = {a1, %2, a1 + ag, % + a2 tas} 
and the basis of p may be taken to have elements i (€¢. + e_«) and (e. — es) for « 


=f 
p,a has basis element 
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hy = i (ee + e-x) for a = {ag, a, + a, %, + a + a3 + a}.  ...(6.5a) 
And the basis of m then can be taken as 
— th, = (ea + e_«) for « = a). ...(6.5b) 
The /;’s give 
Ay = {— %, %, — ag, aa, — (%1 + Oy), — (a3 + ay) 
— (a, + a + a3), — (a, + a3 + a4), — (a, + as), 
— (a, + a, + a + a,)} ...(6.6) 
which gives the set p, 
= {— a1, a, — Og, %q, — (a, + a2), — (H%2 + a3), — (Hy + o%4), 
— (a2 + ag + a4), — (%, + a2 + a3 + o%)}. (6.7) 


The elements of n are obtained by V- e. for « € ps. We see that for this superalgebra 
there will be 2 = 8 parabolic subalgebras. 


(6) B(1,1)—The positive roots of g = B (1,1) are a), a, %, + a, 4, + 2a, and 
2a, — i : 


From Satake diagram we get exp {« (h)} = + 1, fora = 01, %, andexp {a (h)} 
=-—1 fora = «, + «,a, + 2«,, 2%, + 2«,. Therefore k has a basis of ih. for «=a, 
Oo, Cx + Cx, i (Cx — e-x) fora = a,, a. pconsists of basis elements i (e. + e-x), 
(€. — @-x) for « = a; + a, a, + 2x2, 2%, + 2a. In the vector space p, athe maxi- 
mal abelian subalgebra has basis element, 

hi =i ea + e-«) for a = {2a, + 2a}. ...(6.8) 
m is also one-dimensional with 

— ih’, — (€. + €_«) for 4 =. ...(6.9) 
It follows that V = V. with « = 2a, + 2a, and 


—_— — 2 $ 
a. ea Oe, PR TR, PA oy tee a) 
fh | (2%, + Que, 2a. + st [2he, 2] 


and 


Foe { Tea pies (6.10) 


The positive root space A* is obtained from eqn. (4.3) as 


At = {a1, He, & + a2, % + Ze, 2a, + 2x} ...(6.11) 
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which gives set p, using eqns. (4.4) and (6.11) 
P+ = {a,, Oe, &y + He, %, + 2a2, 2a, i 2%}. ..-(6.12) 
The basis elements of n are given by 


$ {€(2014209) rs €_(20%+20,)} =< 3 [2he, 7 2he.); 3 (Ca, 4 €-«;) 


+ 1271/2 g9n N — a, a, + a 


X (Cajtag + €-(ajta9)) + 4 sgn Nasa, N-asstay (Cup p20 

+ €-(a,42%9)); $i (€x, — €_x,) — 4 sgn N-wos0, 409 

X (Cajsag — €-(x,40,)) + $i sgn Noo 1% 1404 N—wosny +09 (a +2% 5 
+ €(a:42%)), $(€x, — €-«,.) —} i sgn Nay,%9 (€a;400 — €(a14%9))5 


$ i (Cas as C_uo 2 5 $ sgn Na, 9%o (Cx 1409 + C(x +09). 
...(6.13) 


This superalgebra contains 2! = 2 parabolic subalgebras. One is the minimal 
parabolic subalgebra given by eqns. (6.8), (6.9) and (6.13) as P: = mM @ a @ nand the 
other being the superalgebra itself. 


(c) B(1,2)—The positive roots of g = B (1, 2) are a1, &, as, a1 + a, & + ag, 
Ge Hb 2us, 2x2 + 2exg, 1 + % + os, 0, + a + 2a, %1 + 2x. + 2as, 2a, + 2a, + 2a. 
From Satake diagram we get 


exp {a (h)} ne ae ] for o= {a1, Xo, Xs, %) a Xo, X, 4 Xe + a3, 2a, + 2a. 


+ 2a3} ...(6.14) 
and 


exp {a (4)} = — 1 fora = {a + a, + 2az, % + 2a, + as}.  ...(6.15) 


k has a basis ih, for « = a, He, X3, € + €5 andi (e, — €_«) for a’s given by eqn. 
(6.14) and p has basis elements (€. — e-«) and i (e. + e_.) for »’s given by (6.15). 


a is of dimension one with 


Ra = {2u, + 2«,} and the element is given by 
hi =i (ex + €-a) for « = {2a, + 2a3} --.(6.16) 


and m is given by 


a7 hit = (ea AG e_a) for wh {a3, 2a, = 2a, + 2a3}. elGch a) 
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The root space A* is given as, using eqn. (4.3), 
At = (—a;, Bo, 3, — (a + do), he + X3, a, + 2u3, 2a, as 2a, @1 + 


+ a3, a1 + ay + 2as + 2a, 2as, 2a, + Do + 2a5}. ...(6.18) 
The subsets p_ and p, are given by (using eqn. (4.4)) 
P— = (a + 203, 2%, + 2a, + 2x5} 
P, = {—a1, Og, &, — (a, + Xo), a2 + ag, 2a, + ag 
ti + % + as, a + a, + 2a, % + 2a, + 2a}. ..-(6.19) 


The elements of n are given by 


° l 
ae {Cras 4205 oe C~ (204205) — $1 (2he, + 2has); aoe (Cx, + e-a,) 


I 
Sgn Naa, {Cx, 42054203 + C- (Hr 42054009)}3 Carpog + C-(w 409)5 
/2 1 1 3 2 ; 
i 


] 
ey (xs = Cx.) + oe sgn Nases (Cas +2x5 — €_(u542%5)); 


i 
AFD S80 Nasxrta3 (Cxgtag — €-(ay405)) + V2 (Castes — €(x9405))5 
(Cu +02 +05 oe C_ (x, 409403); ear teshan, 3% C(x, +r94205)}5 


i 


l 
Vie: (Cx +2%5+2%3 — €- (a, +2%9+208)) oF v2 (sen Na, —(a.,+2%2+2%5)) 
(Gu, F exa,). _..(6.20) 


The minimal parabolic subalgebras is given by eqns. (6.16), (6.17) and (6.20) as p;= m 
@® a @ n. This superalgebra has got two parabolic subalgebras, the minimal parabolic 
subalgebras and the Lie superalgebra itself as the dimension of a is one, 


(d) B (2,2)—The positive roots of g = B (2, 2) are 
1, Lo, Ag, Oy, Bi + XM, to + Ay, as + Hy, Xs + Zag, a1 + % + as, 
he + as + a4, %, + a, + x3 + a, a + ag + 2ay, Zao + 2axsg + 2a4, 
ape ay cde aa rt eg. 1. 2 x, eed Jay) 4-2e,, 
a, + 2%, + 2a, 2aq + 2%, + 2as+ 2a, 
From Satake diagram we have exp {a (4)} = + 1 for 
& = {H1, Oy, Hg, Xy, Xe + Og, & + Hy, %y + Hy, Hy + Hy + Ay, Ay + My 


+ ds + O4, 2a, = 2a + 2a + 2x,}. ...(6.21) 


1074 VEENA SHARMA AND K. C. TRIPATHY 
and 
exp {a (h)} = — 1 for a = {a1 + @ + 2a3 + 204, a, + 2a. + 2a 
A eg, + og + Oy, Hy + Zag + Zag, H+ He + as + 2g, My 
+ a3 + 2a4, 2a. + 2x3 + 2ay}. eet Garo hk; 


The basis of k is given by iha for « = 0}, %, a3, %q, (Ca + €-«), t (Cn — e-«) for « 
given by (6.21). Similarly the basis of p consists of elements (¢. — e-«), / (€, + e_«) 
for « given by (6.22). 


The maximal abelian subalgebra a has one element 

hi =i (ex + e-x) for « = {2a, + 2a3 + 2a}. OL ear 
The subalgebra m is three-dimensional given by 

— ih, = ihe for B = {us, a3 + 2a,g, 2a, + 2a. + 2as + 2a4}. ...(6.24) 


Therefore we can write A* as, using eqns. (6.23) and (6.24), 
At = {—«a,, — G2, — a3, — O, %, + a, — (a + a3), — (as + ay), 
— (a3 + 2a4), a + a% + a3, — (a. + as + a), a +a, + a, 
+ (a4, — aH + 2a + 204), — (a + a3 + 204), — (20, + 2us 
+ 2a4), % + ae + Oy + Zag, a + ae + Zag + 204, ar + 2a, 
+ 2a3 + 2a4, 20, + 2a, + 2as + 204} .-.(6.25) 
Subsequently the two subsets p, and p- can be written as follows: 
P- = {% + Oy, %1 + ae + ag, % + GH. + Hs + a, — (ae + 2x3 + 2a), 
a, + a, + 2a, + 2ag, 
Ps = {—G1, — Gy, —A3, — ay, —(% + as), — (as + a4), — (as + ars) 
— (42 + a3 + %), 1 + Zug + Zug + 2ay, — (ay + 2a%3 + 2a,), 
X ay a + as + 204, 2a, + 2a, + as + 201,}. .- (6.26) 


The subalgebra n is characterized by the algebra spanned by elements V-1¢, for all 


“% € px given by eqn. (6.26) and the elements of n are defined as n = n () g. The mini- 
mal parabolic subalgebra is given by eqns. (6.23), (6.24) and the elements of n. This 
superalgebra also consists of two parabolic subalgebras only, one the minimal parabo- 
lic subalgebra and the other superalgebra itself. 


(e) D (2,1)—The positive roots of g = D (2, 1) are given as «,, a, as, a. + ho, 
Hy b as, 1 + ao + as, 2a, + a, + as, For this superalgebra 
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exp {a ()} = + 1, fora = {%1, %9, %3, a, + ay} ... (6.27) 
and exp {a (4)} = — | fora = {0 + a3, ou + a 4 M3, 2a, + a» + as}. 
...(6.28) 


The basis of k consists of elements ih, for « = Gs %, and (en + e_.), i(e, — €_«) 
for « given by eqn. (6.27). And the basis of p consists of (e,, —e€-«), i (ex. + e_.) for 
% givea by eqn. (6.28). The Subalgebra a is one dimensional and has basis element 


hy =i(e. + e-.) with x = {2a, + a + as) (6.29) 


and m is given by 
th, = ihe for B = {ap, as}. .--(6.30) 
The /,’s give the positive root Space A* as 


A+ = {a, SS rere Xs, (a; ae as), Hy TS AZ, KH) a¢ Xo 9 a3, 2a, a hs + a3}, 
...(6.31) 
The two subsets p- and p, are given by 


P- = {a,, — do}, 
en eur ays cet ae tag, Dart on bah 4..(6.39) 
The elements of n are given by V-l e, fora € P.andn = ng. This superalgebra 


has also two parabolic subalgebras. one being the minimal parabolic subalgebre and 
the other being the superalgebra itself. 


(fee D: (2:2)—The Positive roots of g = D (2,2) are a1, Ho, a, ha, 1 + a, Oy 
a %3, H + Gg, 1 - a, + %3, = 5 as + a, % + XH + ay, 2a2 + a3 + a4, 20 
+ Qu, + a3 + a, a, + % + a3 + a, % + 2a, + a3 +a, The Satake diagrams 
give 


exp {a (A)} = 1 fore = {a1, he, H3, hy, XH, + Xo, He &3, % + ha, 2a, +h 


cea} (6,33) 
and 
exp {a (hA)}} = -— 1 fora = {a1 + a, + as, te + oy + Oy, a, + he + XK, 
20, H 2a, + as + ey, 1 + og + wg + Hy, oy H+ 2%. + %3 + ca}. 
»+(6.34) 


The basis of k consists of ih, for « = H1, Moy Og, Xa, (Cn + Cx), i (Cn — e_g) for « given 
by (6. 33) and therefore the basis of Pp is given by (€, — e_a), i (e, + €-«) for « satis- 
fying (6.34). The basis of a is 


h, =i(e« + ea) for « = 2a, + 2% + a3 + a, oss(0.35) 
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and the basis of m may be taken containing elements 
— ih, = ihe for 8 = {a3, &4, 242 + a3 + a}. ..-(6.36) 


The root space A+ is obtained as 
At = {—a,, —G@, — Gg, — %, — (a; + %), — (a2 + aa) — (% + a4), 
— (2a. + a3 + a4), — (4) + ae + a3), — (a2 + % + a4), — (a 
Tat Oe a4), — (20, + ho + Zag + a4), (a, + a% + 2a, + 2 


as + m%, — (4, + 2a, + Hy + %)}- .. (6.37) 
The two subsets of A+, p, and p_ are given as follows : 
p. = {— (4 + Si). oe a4), — (4, + a + a + «4)} ...(6.38) 


Pig A Sy gy tig ys (a2 + tg), (2%_ + a; + Ga) — (4,4 & 
+ 3), ad (%2 + a3 + a), — (a -- Lo + Oe), =< (2a, + 3 + Os), 
— (% + 2a, + a3 + %,)} .. (6.39) 


~ 


The elements of n are given by V-!e, for « € pt and n= nf) g. This superal- 


gebra will have two parabolic subalgebras, one the minimal parabolic subalgebra and 
the other the superalgebra itself. 


6. CONCLUSION 


We have presented here how to obtain the parabolic and minimal parabolic sub- 
algebras for some Lie superalgebras of physical interest. Schmidt construction, it is 


hoped, will yield the various representations. We will report more on specific cases in 
a forthcoming communication. 
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SOME PROPERTIES OF THE SPHERES IN METRIC SPACES 


THOMAS KIVENTIDIS 


Department of Mathematics, University of Thessaloniki, Greece 
(Received 5 October 1987) 


In a metric space (X, d) if for any distinct points x, y & X there exist seque- 
nees (s,) > y, (s\) + x, ny, s Fxyn € Nwith 


max {d (x, Sn), d(y, 5n)} < d(x, y) 


and 
max {d(x, s'),d(v, 8° )} < d(x, y) 


then the closure of a open sphere is equal to the closed sphere. 


Furthermore, if we suppose that for any distinct points x, y € X there 
exists z € X such that B(x, d(x, y)) 1 B (z, d(z, y)) = ¢ or the metric space 
is externally convex”, the interior ofa closed sphere is equal to the open 
sphere. 


§ 1. Let (X,d) be a metric space, B(x,r) = {y © X¥: d(x, y) <r} the open 
sphere, B (x, r) = {y © X: d (x, v) Sr} the closed sphere, clA the closure and int4 
the interior of asubset A C YX. 


It is well-known® that in a normed space the following properties hold : 
Property A: for every x © X,r > 0,cl B(x, r) = B(x, r); 
Property B: int B (x, r) = B (x, r). 


The properties (A) and (8) do not hold in all metric spaces, for example, if (X, d@) 
is a metric space with d (x, y) = 1 for x # y and d (x, y) = 0 for x = y, then subsists 
B(x, 1) = {x}, B(x, 1) = X and Cl B(x, 1) = {x} 4 B(x, 1) = X, int B (x, 1) 
= X + {x} = B(x, 1). 


§ 2. Given any two points «, B, « « $8 of the metric space (X, d) we put E (a, 6) 
= {z © X: max {d (a, z), d (8, z)} <d (a, B)} U {«, B} and we denote by E' (a, B) 
the derived set of E («, 8), i. e. the set of accmulation points of E («, 8). 


Lemma |'—For every x € X andr > 0, we have cl B(x, r) = B(x, r) if and 
only if for any «4,B € ¥,a #8, E (a, B) C E’ (a, B). 
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Proposition 1—If in a metric space (X, d) the following condition holds : for any 
distinct points x, » © X there exist sequence (s,) > y, (5, ) > %, Sn FY 8, FX 
¥n€E N with 

max {d (x, 5,), d (Y, 5,)} < d (x, Y) 


and 
max {d (x, s,), d(y, s,)} < d (x,y) 


then for every x € X,r > 0, cl B(x, r) = B(x, r). 
Proor : According to the lemma we need to show that any distinct points 
a,8 € X, E(a, 8) C E’ (a, 8). 
Letz € E(a, B),z~4,z 4 (clearly «, B € E (a, B) by hypothesis). 
For the points «, z € X there exists a sequence (s,) > z with 
max {d (a, s,), d (z, S,)} < d (a, z). 
Observe that for 0 < « < min {d (a, B) — d(a, z), d (a, B), — d (f, z)} exists 
an integer ™) = mo (e) € Nsuch that: ¥ n => nm, 5s, € B(z, «) and 
d (s,, 8) Sd (s,, z) + d (z, 8) < « + d(z, B) < d(a, 8) 
d (sn, %) S d(s,, z) + d(z, «) << ¢ + d(z, a) < d(a, 8). 
Hence, 
s, © E (a, 8), form =m. Thus z € E’ (a, 8). 


Remarks 1: (i) If the metric space (X, d) is convex in the sense of Menger? (i.e. 
for any distinct points x, y € X there exists at least another point s © X such that 
d(x, s) + d(s, y) = d(x, y)) and complete, then the hypothesis of Proposition | 
holds (see Theorem 14.1 in Blumenthal’). 


Generally, if the distance d is convex (i.e. for any distinct points x, y € X there 


exists a point z € X such thatd(x,z)= d(y,z) = aay the hypothesis of 


Proposition 1 holds. We note that the convexity of the distance ddoes non imply 
that the space X is complete, for example, if ¥ = (0, 1) with d(x, y) = | x — FAS 


(ii) Let X be the space of the irrational numbers of the open interval (0, 4) and 
the rational numbers 1, 3 with the distance d (x,y) =|x-—y|. 


It is obvious that the distance dis not convex in X. But, it is easy to prove that 


the hypothesis of Proposition! holds and the space Xis neither complete nor 
connected. 
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(iii) One verifies easily that every metric subspace (Y, d), Y C X of a metric 


space (X, d) which has the property (A), with Y an open or everywhere dense subset of 
X, has also property (A). 


§ 3 Let X be the closed unit sphere in the plane R® with the Euclidean distance 
d,. Evidently, for every x € X,r > 0, B(x, r) = cl B(x, r). In particular, for x = 0 
(the origin), r= 1 we have B (0,1) = X¥ = cl B (0,1) whereas int B (0, 1) 
= X¥ + B(0, 1). 


We can verify that in the space X¥ the hypothesis of Proposition 1 holds (the 
distance d» is convex in X), 


Proposition 2—In a metric space (X, d), which satisfies the hypothesis of Propo- 
sition 1, if furthermore for any distinct points x,y € X there exists another z € X 
such that 


B(x, d(x, y)) O B(z, d (zy) = ¢ 
then for every x © X,r > 0 int B(x, r) = B(x, r). 


Proor: Let B(x, r) be any closed sphere. We need to show that no point 
y € X with d (x, y) = ris an interior point of the closed sphere B (x, r). 


In fact, by hypothesis there exists z € X such that 
B (x, d (x, y)) O B(z, 4 (z, y)) = ¢. 


Since the metric space (X, d) satisfies the hypothesis of Proposition | there exists 
a sequence (s,) > y with 


max {d (z, S,), 4 (y, S,)} < d (z, y). 
We have d (s,, z) < (z, ») thus s, € B(z, d(z, y)) ands, & B (x,d (x, y)). 


Then s, & B(x, d (x, y)) because, if d (x, s,) = d(x, y) (for some s,) then since 
s, © B(z,d (z,y)) there exists a neighbourhood B (s,,€) contained in B (z,d (z,y)) and for 


the points x, s,, again by hypothesis of Proposition 1, there exists a sequence (s\ )—> Sp 


with 


max {d (x, 5’), d (Sn, 5, )} < d(x, Sn) 


and we would have s, € B(x, d(x, y) NO B, d(z, y). 


But this contradicts the hypothesis. 


Hence, in every neighbourhood of the y there exist points which do not belong 
to B(x, r), i.e. y & int B(x, r). 
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Definition 2—A metric space (X, d) is externally convex provided it contains for 
each pair of distinct points x,y € X at least one point z € X such that d(x, y) 


+ d(y, z) = d (x, 2). 

Proposition 3—If a metric space (X, d) satisfies the hypothesis of the Proposi- 
tion | and furthermore, is externally convex then for every x € X,r > 0, int B(x, in 
= B(x, r). 

Proor : We will show that no point y € X with d (x, y) = ris an interior point 
of the closed sphere B (x, r). 

Since the metric space is externally convex there exists one point z € X such 
that d(x, y) + d (y, z) = d (x, z). 

By the hypothesis of Proposition |, there exists a sequence (s,) > y with 

max {d (y, $,)}, d(z, 5,)} < d(y, z), ¥n EN. 
We need now to show that s, & B(x, d(x, y)). 
In fact, if s, € B(x, r) then we would have 

d (x, Sn) S d(x, y) > d (x, 5,)+d (S,, z) < d(x, y) + d(y, z) = d(x, z). 
But this contradicts the triangle inequality. 


Remarks 2: (i) Consider as metric space (X, d) the positive angle of the plane 
R’ Oxy including both the positive semiaxis Ox, Oy, with the Euclidean distance. 


Observe that in XY we have int B (x, r) = B(x, r) whereas the conditions of the 
Propositions 2 and 3 except for the hypothesis of the Proposition 1, do not hold for 
the points x € Ox, y € Oy. 


Consequently, the conditions of Propositions 2 and 3 are only sufficients. 


(ii) One verifies easily that each metric subspace (Y,d), Y C X of the metric space 
(X, d) which has the property (B) by virtue of the conditions of Proposition 2 or 3, 
with Y an open or everywhere dense subset of Y, has also the Property (B). 


(iii) Consider as metric space (X, d) the open unit sphere of the plane R? with 
the Euclidean distance. Then the conditions of Propositions 2 and 3 hold, whereas if 


the space XY is the closed unit sphere these do not hold, except the hypothesis of 
Proposition 1. 
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EFFECT OF PULSED LASER ON HUMAN SKIN 


D. RAMA MurTHY AND A. V. MANOHARA SARMA 


Department of Mathematics, University College of Science, Osmania University 
Hyderabad 500007 


(Received 22 April 1987) 


Taking skin asa semi-infinite homogeneous solid with distributed source 
along a finite depth from the surface, the temperature and thermal stress dis- 
tributions are obtained, using modified hyperbolic heat conduction equation. 
The effect of pulsed laser on human skin is considered. Parameters of phy- 
sical interest are plotted. 


1. INTRODUCTION 


Lasers find wide use in numerous fields of Science and Engineering. The laser 
produces radiation with a highly regular light field, outstanding in its high coherence, 
monochromaticity and directivity. Since laser beams are of high power and can be ex- 
ploited to produce a targetted effect on material. The applications of lasers include 
welding, cutting hole burning, isotope separation and medical diogonising etc. Further 
more, it finds applications in data transmission and processing, measurements and 
quality control. 


The tolerance of the human body to microwaves has been studied by Hendler 
et al). The skin of the human body responds differently for different wavelengths of 
electromagnetic radiation and as such the spectral characteristics of the laser emission 
will determine the tolerance. 


In this paper, the effect of pulsed laser on human skin, using modified heat con- 
duction equation is considered. The analogous temperature and stress distributions, 
which includes the effect of finite speed of heat propagation is determined using Laplace 
transform. The tempearature distribution is computed and exhibited graphically for a 
ruby laser. The results due to Majumdar’ can be obtained as a particular case. 


2. FORMULATION AND SOLUTION OF THE PROBLEM 


The simplest physical model closely resembling normal skin is feasible to mathe- 
matical analysis as a translucent and semi-infinite solid. Following the principle of 
conservation of energy, shorter the pulse duration, the smaller will be the depth of 
penetration by the radiation, for a given energy input. It is also reasonable to assume 
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that the thermal inertia of the tissue will prevent any significant change in its thermal 
and optical constants during the shorter exposure period. 


Consider the skin as a semi-infinite homogeneous solid with distributed source 
along the depth from the surface. Initially, the temperature is considcred to be uni- 
form. The differential equation of heat conduction in this case is 


1 &T ter 2T q 
Gath eo eee pC, ---(1) 








where C is the velocity of propagation of heat, T the temperature distribution, A is the 
thermal diffusivity, g the rate of energy absorption per unit volume, P the density and 
C, the specific heat. Equation (1) is obtained from modified heat conduction equation 
(see Chester*, Baumeister and Hamill*) by considering the finiteness of the heat pro- 
pagation velocity. The initial and boundary conditions are 


at 


oT 
t = 0, T(x, t) = 7, and eee 0. TA C4 


Here 7, is the constant temperature. 


The regularity boundary condition is T (eo, t) = 0 and the effect of penetration gives 
at 


6T 
x=0 ra = Hf (] —r) ee kh 


where k is the thermal conductivity, H is the constant intensity of the incident radia- 
tion and r is the spectral reflectance of the skin for a particular wavelength (A,). The 
right side of equation (3) determines the rate at which energy is absorbed per unit area 
of the surface and the rate of energy absorbation (q) is 


q = «H e-*, (4) 
In eqn. (4), « is the absorption coefficient. 


Employing the nondimensional variables 


C =< 
Zi. seek Se Plena lk 
oh om) Rie eee 7 





the heat conduction equation (1), the initial and boundary conditions (2) and (3) can 
be written as 


89 00 629 
5p? + 2 Fy Wier Pe 4:9 me te 


at 
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Pay fea) ee F, = 0, 0(0%,B) = 0 6) 
at 
Meee a dO (1 — 7) 
z= 0, So Ae) 
Here 


Q = 4h? «H/C? Tok, s = 2h/C. 
Applying Laplace transform to (5), (6) and (7), one obtains 


Gl gery eee eee 
sh pvp? + 2p PV p*+2p[(p+1)?—(1+s")| 


pee 
= Pl[(pt+1)?. — (1+s?)] ..(8) 





6= 


where 


@ is the temperature in the transformed domain and p is the transform parameter. 


The inverse laplace transform of eqn. (8), from Bateman’, is 


B : B 
@=A J gi (u) du — BS 8 (u) 81 (B — u) au 





B 
ie ae | e~“ sin h (Ru) du Fe be 
0 


where 


Q (lar) Ss 
oe sh Se 


g, (2) = 0, if B < 2 
= eF 1, (6 — 2), ifB > z 


_ e® sin h (RB) 
82 (8) a R : 


Here, J, is modified Bessel function of zero order and R? = (1 + 5s”). 


The temperature distribution given in eqn. (9) is evaluated numerically. 


For the determination of normal thermal stress ox:, consider the following dyna- 
mical equations of themoelasticity. 
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da 02u 
me FG ais ... (10 
Ox p or? (10) 
ou 
Gxx = (A + 2n) 7 - — BT ad OW 
where u is the displacement, A and » are familiar Lame’s constants and B, = (3A + 
2u) A. 
From (10) and (11), we get 
Ou 1 Ou oT 
Ox vets DER ee ast 
Cy 


Here C; = V (A + 2u)/P is the velocity of propagation of longitudinal wave and 
m = B,/(A + 2u) 


Introducting the potential of thermoelastic strain ¢, where 





od 
uti=> ax: 3 ati) 
Using 
C é: 
U = aR M and z = >A x We get from eqn. (13), 
_ 9" 
U= Pe ...(14) 
where 
CT 
eT 


From eqns. (12), (13) and (14), we get 


0*g* a7h* 
agt —  —apr— = m8. BA Gm)" 





Here 
at = CUCrs 


Applying laplace transform to (15), we obtain 


$* = ST eo ee: e-Vo's0p 3 | 
ae 2 — 
P’ (p+2E) [(p+1)?? — R* ap Vp? + 2p 





. Sues om aie ea tack dr 
D* (p + 2E) ap qe 
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7 OS | e-tpz ae 
pip —R)(e— bh) ap ss ip 





...(16) 
Here, 


Pe 
a 





l1— @ 


Following the method of superposition as given in Nowacki®, the normal stress o:2 is 








= c= = he a b* a2 * 
Ger = Ose + Or2 = P ( a aa ) 17) 
Here, * is determined using boundary condition 
oz; = 0 at z = O and the equation 
2 ),* 2 |)* 
Aah Meileae. pce dinar t (18) 


az2 9 OB? 


Using Laplace transform, from eqn (17) and (18), we obtain 





te — [ mBE mee l ] 
y: p(p + 2E)((p +1)? —- RF] Lap Sp?+2p 
pee abies feb al | 
p?(p + 2E) L ap Jp? + 2p 





is a’p[(p + 1)? = R*] By | ap me = il eam meee ED 
and 
Ger = P pt[o* + 4]. ...(20) 
Now, from eqn. (16), (19) and (20), we obtain 





= mBE ee Ap aa 
On = P Peace iee Rr [e é 
(p + 2E) ((p +1)* —R'] vp* +2p 
mAE e- 2 eV p*+op z | 
(p+2E) L /p* + 2p Jp? + 2p 
pes SOE ceca OE ae QI 
+ a pi(p + 1)? — R*] (p® — 8%) s s a 


From Bateman’, the inverse transform of eqn. (21) is 


o,; = P [mB (F, — Fy) — mA(F, — Fu) + Q (F, — Fo) ...(22) 
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67459 
4248-1 


2749-2 


(%) 


mnNcdrnmozme 


283 .74 


283.71 


283 .60 





ToRME ECR) 


Fic. 1. Temperature Vs time. 


where 


B B 
Bom BLD OO — a) fib ity Fin Eas ee ae 


B B 
F; =E J e~ shy U(y —. az) fy (8 a y) dy, ha Es extky f, (B — y) dy, 
1 ; 7 
Fy = =| U(y — az) f, (B = y) dy, Fz = | fe (— y) dy 
a 
0 0 
B B 
h® =|LOKC-Yo,L@ = —| 2" 
J 7 y(B — y) 
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p=10 


4747-70 


4746-95 
3248-91 


(%) 


mrpcaremvem 


ow 
iS 
> 
@ 
oO 
oO 


1250.39 


1249-84 


ee Ca ee UEEEEEEEEEn ienenendtl 


8 41 15 
DISTANCE (2) 
Fic. 2. Temperature Vs Distance. 





B 
fr (8) = be ® sinh (28), f. (®) = | fe 0) e8® dy, 


B 


i < 6% WEP) UC—2, 40 = 7 | sab ov 


fs (B — y) dy 


in the above expressions U (y — az), U(B— z) are unit step functions. 


3. RESULTS AND DisCUSSION 


In order to interpret, the temperature variation, the initial temperature of the 
skin is 33°C, PC, = 1, k = 107-*c. g. s. units are taken. The values of optical con- 
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stants depend on the nature of radiation employed. In addition, for computation, we 
take ruby laser of wavelength A; = 6943°4 with the following values for the parame- 
ters involved, H = 0.1, r = 40 and « = 20. 


The temperature distribution given in eqn. (9) is evaluated numerically for the 
parameters given above. The variation of temperature distribution at z = 7, z = 10 
and z = 15, respectively are presented in Fig. 1, where as in Fig. 2, it is given for 
6 = 3,8 = 7 and B = 10. 


In Fig. 1, it is seen that the temperature is linear one and increases with time, in 
general, and as the rate of penetration increases, the temperature will also increase. 
This phenomena is practicably feasible. In Fig. 2, temperature as a function of dis- 
tance for different values of time (8) is plotted. As the time increases the temperature 
will also increase and however for higher values of distance, the temperature keeps on 
decreasing. As the time increases, the slope of curves decreases. 


These results included the effect of finite speed of heat propagation. We can ob- 
serve the inherent wave nature of the heat transfer process and the temperature at any 
depth is always Jess than that of surface temperature during the exposure period. From 
eqn. (9), the results due to Majumdar? can be obtained as a particular case. 
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In this paper, an arbitrary function f(x) defined on the interval (0, a) is 
expressed as an expansion in Dini’s series of the orthogonal family 
{Fy (em x)} of modified Bessel functions, where Gy (x) = x’ Jy (x) and pm 
denotes the mth postive root of the equationx J", (ax) +h Gy (ax) = 9. 


Next, the convergence theorem is rigorously established. The Dini’s series 
suggest to consider a varaint of the finite Hankel transformation, which will 
be called the finite Hankel-Schwartz transformation of the second kind. This 
transformation is used in solving some partial differential equations which 
cannot be directly treated by applying the corresponding finite Hankel trans- 
formation. Finally, we remark that the initial term of our expansion depends 
only on the parameter h, whereas the classical Dini’s expansion depends on 
h + v. 


1. INTRODUCTION 


Schwartz? investigated the following modified Hankel transformation 
oo 
F(y) = 5 F (xy) f (*) dm (») mA ict) 
0 
Shae dnianl2 oct lr x7. dx and SF (x)iei20T Gor 1) x-¥ Jy (x), Jv (x) 
being the Bessel function of the first kind of order v. 


This transformation has been studied in spaces of distributions by several authors 
and Lee’ calls it Hankel-Schwartz integral transformation. 


To consider the expansion of an arbitrary function f (x) defined in the interval 
(0, a) asa Fourier-Bessel series, i.e., as a series of the type 


f= E a, Fo(in) (1.2) 


where Fv (x) = x-Y Jv (2) vee = 4 and j/, denote the positive zeros of the functions 
Fy (ax), i.€., 


Fu (jn) = 0. eg 
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Méndez? introduced the corresponding finite transformation through the equation 
yy [f (x)) = Fv (1) = J x?%*2 BZ (J, x) f (x) dx ...(1.4) 
0 


which is called Hankel-Schwartz transformation of the first kind of order v. Its 
inversion theorem is stated as: 


Theorem 1—Let f (t) be a function defined in (0, 1) and assumed to be absolutely 
summable over the same interval. Let v > — } and 


1 
a, = 2 | wt (4,1) f(D dt, = 1,2, 

A} ¥*, Ch) 0 
If f(t) is of bounded variation in (a, b),0 <a<b< 1, and if x € (a, b), then the 
series (1.2) converges to 


a[f(x + 0) + f(x — 0)]. 


In this paper we show how the modified Dini series expansion of an arbitrary 
function f (x) leads naturally to the finite Hankel-Schwartz integral transformation of 
the second kind. The inverson theorem of this new transformation is rigorously 
established by studying the convergence of the series expansion. The operational 
calculus generated is used in the solution of several problems in Mathematical Physics. 


Recall that the form of the Dini series is determined by the nature of the zeros 
of the equation 


ZF) (2) +h F(z) = 0. 
We emphasize that the first term of the expansions depend on the parameter h, and 
not on A + vy, as it happens in the classical theory (cf. Watson?*, p. 597). 


Another interesting feature of this transformation is its usefulness in the solution 
of partial differential equations which cannot be tackled by applying the corresponding 
finite Hankel transformation when v # I (cf. Colombo‘, p. 82). 


Finally, we will denote in the sequel the partial sum of the series (1.2) by 


S, (x) = = am Fy im x). (1.5) 
Setting 
P, (x, t) = 2 Fs (Im X) Fv (jen t) .--(1.6) 


ont ie oy me 
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we have 


Ss.) = j 124) P(x, t) f(t) dt. (1.7) 


2. PRELIMINARY RESULTS 


Let L denote the differential operator x~*¥-! é xures <. We begin by consi- 


x 
dering the following Sturm-Liouville problem (cf. Sneddon", p. 440) 
(L+A*°)y=0,0<a<b ZAZA) 
My (a) = a, y (a) + a, y’ (a) = 0, Ny (6) = b, y (6) + by y’ (6) = 0 
mee) 


where a,, a,, 6, and b, represent prescribed constants. 
The general solution of eqn. (2.1) is 


y = $(x,A) = A (A) Fy (Ax) + BOA) WAX) 5 (2:3) 


where Fy (x) = x~-¥ Jy (x) and M (x) = x” Yy (*), Jv (x) being the Bessel function 
of the first kind of order v and Y, (x) the one of second kind. 


Let y = ¢, (x) be the eigenfunctions of the problem (2.1)—(2.2) which corres- 
pond to the nonzero eigenvalues A,. We have the general orthogonality condition 


fore | 2" & eC x82 Oo) 


+ 29 44(x) #% GO} |, ifm =n 


| 
q 0 sl ieee, 
ie 


f x41 by (x) bm (x) dx = 4 
a 


Then, we deduce from (2.3) that the solution of the particular problem 
(L+ *%)¢(x)=0, 0<x<a Te), 
N¢ (a) = ¢' (a) + he (2) = 0, A > O 
is 
gn (x) = Fv (P, x) ...(2.6) 
where P,, P,, ..denote the positive zeros arranged in ascending order of magnitude of 
the transcendental equation 


z F' (az) +h Fy (az) = 0 ...(2.7) 


i.e., [cf. Gray et al.®, p. 16, eqn. (24)], 
—az? Fy, (az) +hFv (az) = 0. 
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The above orthogonality condition (2.4) now becomes 


+y 
( <x (ah? + aP® —2vh) F2 (Ppa), if m=n 
| x2vt1 Zy (P,,x) Fo (Pax) dx = Pn 
F L 0 , if mn. mes) 


Given an arbitrary function f (x) defined in the interval (0, a), (2.8) allows one to 
formally express this function as a Dini expansion, as follows 


f(x) = E by Fi (Pm x) (2.9) 


where 


2p? fe 
bm = |x (0,3) f (2) dx 
qa?vti Ar Sg (Pma) (ah? 4+ ap’ = 2vh) 0 (2 10) 


m = 1, 2, ..., Pm being the positive roots of eqn. (2.7). 

Note that we can extend the Dixon theorem (cf. Watson’’, p. 480) to the zeros 
of the function (2.7). Indeed, it can be proved that the zeros of the equations 
Ax F\ (x) + BF (x) = Oand Cx F¥* (x) + DFy ix) = 0 are interlaced, what- 


ever the real numbers A, B, C and D, provided they are such that AD + BC. Hence, 
the roots of eqns. (1.3) and (2.7) also are interlaced. 


3. THe Mopiriep Dint EXxPANSION—THE CONVERGENCE THEOREM 


In this section it will be assumed that a = 1 for the sake of simplicity. 

As it occurs in the classical theory of Dini expansions (cf. Watson’’, p. 597), we 
have to add an initial term to the series (2.9). In fact, the form of Dini expansion is 
based on the zeros of the function (2.7) and these depend upon the values of the para- 
meter h. Thus, the expansion (2.9) only corresponds to the case h > 0, 


When fA = O it can easily be seen that the equation (2.7) has a zero at the origin. 
On the other hand, when / < 0 this function has two purely imaginary zeros. 


Let v be a real number such that vy > — 4. We write the modified Dini expan- 
sion of f (x), as follows 


f(x) = bo +E bm Fy (Pmx) (3.1) 


where by denote the initial term which must be inserted in (2.9) asa consequence of 
the existence of these new roots. 
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If h > Othe initial term 6, = 0 and (3.1) coincides with (2.9), But when 


h = 0, taking into account that 


1 
J > ahh Fy (Pm x) dx = Fr4i (P,,) = 0, m= Wr Pee 


[cf. Gray et al.®, p. 16, eqn. (25)] and (2.7), we get 


ee tfJe 429) j x41 f(x) dx. (3.2) 


Finally, if -- Pp i denote the imaginary zeros of (2.7) when h < 0, from (3.1) 
and (2.8) we infer that 


2p; } 
b 5 = —<—  ———————_ | x"*? Fy (iP x) f(x) dx. ...(3.3) 
(P; + 2vh — h?) a (Po 7) 0 
Now, consider the function 
2w Fy (xw) Fy (tw) 


pee ay ew EW) <r iW )S ...(3.4) 
Fy (w) (wF' (w) + h Fy (w)} 


whose poles are the zeros j,, Js)... of YF, (z) and the zeros Pi, Ps ... of z F< (z) 
+h F(z). 
The residues of this function at the first poles are 
2 Fv (in *) Fv Cin) 
Jen SF s0¢ Cin) 
If h > O the residues at the poles P:, Po, ... are 


2p* Fy (Pn x) Fv (Pn t) 





(p> — 2vh + h*) FY (Pn) 


When h = 0 we must moreover consider the residue at the origin, whose value 
is — 4(v + 1). 
When ht < 0 the residues at + iP, are both equal to 


202 Fv (Po xi) Fv (Po ti) 





(pP? + 2vh — h*) FZ} (Poi) 
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By denoting the partial sum of the series (3.1) 


a, (x) = by + = bm Fv (Pm x) ...(3.5) 
and 
ts 2p*. Fy (P,, x) fv (Pm t) 
P, (x, t; h) = Ag (x,t) + —_ rata) 
mai (P*,— 2vh + h?) F? (P,,) 
where 
{ 0, ifh > 0 
coin SO) 
Ay (x, t) = a al (39) 
ate : 2p* Fy (Po xi) Fy (Po ti) 
| ———_____________ ifh<0 
| (5 + 2h — ht) FH} (Po i) 
we can express (3.5) as 
0, (x) = f t+! P, (x, th) f(t) dt. 3.8) 
0 


Now choose D,, such that itis not equal to any of the number jm and p, < D, 


< Pns, and let jw be the greatest of the numbers Jm Which does not exceed D, (ef. 
Watson’?’, p. 598). 


The following expression 


S, (x, tr hh) = 2Fv (jm X) Fv (jm t) 


= Ao (x, t) 
mal Ie Sis (jm) 


n2p? Hy (Pm x) FH. (P,, t) 
SS 9) 


mai (P+ h? — 2vh) F? (Pm) 


where A, (x, t) is given by (3.7), permits to connect the partial sums of the modified 


series of Fourier-Bessel (1.2) and Dini (3.1). Clearly, from (1.5), (1.6), (1.7), (3.5), 
(3.8) and (3.9) it can be deduced that 


1 N : 
pe" See BAS (dt — Ba Sy ual by ae beeen ey 
pees m=1 


= Sw (x) — a, (x). ...(3.10) 
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From Cauchy’s theory of residues we find the following integral representation 
of (3.9) 





D+01 
5, 4 A) = = | pen 2s FAA) Fy MEW) a yyy 3,11) 
pcos Fv (w) {w F, (w) + h Fv (w)} 

Since 

t 

§ wt?et) F, (tw) dt = tt? Fy, (tw) 

0 
[cf. Gray et al.*, p. 16, eqn. (25)], it can be inferred from (3.11) that 

t Me acs { 

2ut 
| 141, (x, tA) dt = 4 ~ ; Ww Fy (xw) Fos (0) gy 
° D =o t Fy (w) {w Fi) (w) + h Fy (w)} 
ec 12) 
As an immediate consequence of (3.11) and (3.12) we have 
eee se Teme) 

os (x, t, h) | < ead ately (2—x—- t) oa | ) 
and 

. t +1 /2 1 

: : Tail fcboeel Ie <i akc 
iG Sa a8) a, Sons (=) =D (3.14) 





0 
where c3 and c, are constants independent of 7, x and ¢. 


Next, it can be proved with an argument similar to the one used in Watson™ 
(p. 599) that if f (t) is absolutely summable in the interval (a, b),0 < a < 6b < 1, then 


i t2¥+1 § (x, t; h) f(t) dt > 0, as n > © Hi ERE), 


provided 0 < x < l. 


Theorem 2—Let f (t) be a function defined and absolutely summable in the 
interval (0, 1). If f(t) is of bounded variation in (a, b) whereO0 <a < b < 1, then 
the series (3.1) converges to the sum A[f(x + 0) + f(x - 0)] at all points x such 
thta+tA<x«qb-—A4, A > 0 being arbitrarily small. 


oo 
; : ies Sa, Fv (j, x) converges to the 
Proor : By virture of Theorem 1 the seri iow v Vj 


sum $ [f(x + 0) +f(* — 9]. Our assertion follows directly from (3.10) and (3.15) 
to pass to the limit as n > °°. 
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Remark | : Note that the initial term b, of our expansion (3.1) only depends on 
the value of h, whereas this term depends on A+ vy in the classical theory 


(cf. Watson’, p. 598). Moreover, the roots P, of the equation z F u(z) + h Fy (2=0 


are not equal to the roots A,’ s of z J; (z) + hJy (z) = (z) = 0. 


4. THe Finrrs HANKEL—SCHWARTZ INTEGRAL TRANSFORMATION OF 
THE SECOND KIND—APPLICATIONS 


According to (2.9) and (2.10), we define the finite Hankel-Schwartz integral 
transformation of the second kind of order v > — } by the equation 


sMoyy [Lf (x)] = Fa,v (n) ea j “ah Fy (P,, Sy dx ...(4.1) 


whose kernel is the modified Bessel function (2.6) and where P, denote the roots of 
eqn. (2.7). 


The corresponding inversion formula is 
2 ee PS Fay (nt) Fv (P, x) 


shy, [Fav (2)] = f (x) = by + goetl NESE ee 
n=1 (ah? = ap* ss 2vh) be i (P, a) 


...(4.2) 


Theorem 2 not only guarantees existence of (4.1) but also ensures that inversion 
formula (4.2) holds. 


If we assume that the f € C? (0, a), f’ (a) + hf(a) = Oandh> 0, we obtain 
the main operational formula of this transformation, i. e., 


” 2v ’ 2 
sMoyv | f (x) + = 5 (x) | = a. Pe sho,v iy (x)] .-.(4.3) 





whatever the values of f (0), and f’ (0), provided they are finite. 
If f’ (a) + hf (a) 40 andh > 0, we get 


rat eae 
shes Ge eae (x) ] = at F, (P, @) (7 (@) + h f(a) 





— Ph shay Lf (x)]. -+( +4) 
Remark 2: Recall that the functiony = GF, (x) isa solution of the equation 


», 1+2y 
Ll= ee ae = iplicati 
ySy’+ nee: +y=0. Its multiplication by x?’ has Only repercussions on 


the sign of the parameter y, that is, the function 


Y= FP) Sates) 
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is a solution of the equation 


Teest2¥ 
x 





Le y=y’ + y +y=0. 


Consequently, the solution of Sturm-Liouville problem 


(L* + A*) ¢ (x) = 0 


N* $ (a) = ¢' (a) + (»- =) 8 (a) =0 


o* (x) = FP P(x) (4.5) 


where P,, P2,... denote the positive zeros of the equation (2.7). The solutions (4.5) 
form a orthogonal system on the interval (0, a) with respect to the weight function 
x1-2¥. Proceeding as before, we can now introduce the integral transform 


Hts FCO = FR = [x FF (6, FC) ax (4.6) 


whose inversion formula, in the case h > 0, is given by 





o0 2p? FX, (n) F* (P, x) 
AS FE, =f) = > | 
SS at (ah? + ap? — 2vh) FY" (Pn a) 


(4.7) 


A similar result to that proven in Theorem 2 can be stated in relation with the 
convergence of the series (4.7), whenever v = = *, 


The main operational rule of the transform (4 6) is 


at [7 + Are) = - 0h A . (48) 





provided that f’ (a) + (1 aes \r@ Spend We 
In the sequel we shall give a few examples to illustrate the use of the above 
transformations in solving some important problems. 


(a) Let v be any real number. We wish to find the solution of the equation 


C’u 2v+ 1 Ou 1 Ou 

es vr Sas = O<r<ak>0 (4.9 
are r or pest ral acca ae ) (4.9) 
satisfying the initial condition 


u(r, 0)=f(r) ©Srs4) 


1098 J. M. MBNDEZ 
and the boundary conditions 


ou (Gs) + h(a, 1) = 0, for every v = 0, 


r 
or 
vial BD, + (1 -2) u (a, t) = 0, for every v => 0. 
or a 
By virtue of (4.3) and (4.8), we convert formally (4.9) into 
a 2 ) v(t) = 0 
(< + ket) tw (0) = 0, 
where 
shoyv [t ir, u)] » Vv = 0 
cede, ie u(r, Ov 0. 
Hence, 
—keir 
Un (t) = Fa,v (n) € ... (4.10) 
where 


sMoyy [f(r)] ree! = 0 
eae a Lf] ¥ <0. 


By applying the inversion |formulas (4.2) and (4.7) to (4.10), we get the required 
solution 











f 7 a p> F,,y (”) Fy (P, r) e~kont 
: ati r eats * ea W 
| n=1 (ah? + apy —2vh) + Mr (P, a) 
u(r, t) = 4 
; 0o p2 F,,y (n) F*. (p, r) enka t sat eLe 
| qae’ti [Sa te pe a ¥ 20: 
t nai (ah° + ap) + 2vh) F*? (p, a) 


Note that when v = 0 the prablem (4.9) reduces to the one considered by 
Sneddon” (eqns. 8-4-20, 33, 34) on the diffusion equation, since in this case Fy (x) 
= Jo (x) and p, = E, are the roots of x J, (ax) +h Jo (ax) = 0, Then, both of 
formulas in (4.1 1) yield the sum solution and this coincides with the one achieved in 
the reference mentioned obove. 


A procedure similar to the one used by Churchill® (p. 191), allows one to esta- 
blish (4.11) as a rigorous solution of our problem. 
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Remark 3: Note that the equation (4.9) cannot be solved directly by means of 
the finite Hankel transformation, except when v = 0. Nevertheless, the simultaneous 
application of the finite Hankel-Schwartz transformations (4.1) and (4.6) provides a 
simple method to solve immediately the problem (a), no matter what the real value of 
v may be. 


(b) Many partial differential equations involving the n-dimensional laplacian 
operator can also be solved by using the transformation (4.1). Indeed, the n-dimen- 
sional potential equation is 











2 2 2 2 

Smee Po OS ES .. (4.12) 
2 2 9 Z 

Ox, Ox, OXi 55 


where u = U (X1, Xo, «225 Xn-1» Z). If we seek solntions which only depend on 
r=x?+(x2+...4 x°2_,)'!* and z, (4.12) reduces to the from (cf. Sneddon’, p. 342) 


=. 


?u n— 2 ou Ou 
re 


ar? ror az? 





=r, (4.13) 


We find now the solution of (4.13) that satisfies the conditions 


UD + hu (az) =0(2 > 0,h > 0) 


u(r, 0) = f(r) ..(4.14) 
u(r, z) > 0, as z > ©, 


by directly applying to (4.13) the finite Hankel-Schwartz transformation of the second 
kind of order v = (n — 3)/2. Now denote U, (z) = sh2,v [u(r, z)]. From (4.3) we 
see that U,, (z) satisfies the equation 


a U, (z 
— P, U, (2) + ah = 0 


whose solution is, in view of conditions (4.14), 
U,, (z) = F (n) e- Pn? 
where 
F (n) = shov(f(r))- 
Again making use of (4.2), the formal solution of the problem posed by equations 


(4.13)—(4.14) is 


cop? e-tnt Fv (Py r) F (n) 





ip e im oe ...(4.15) 
aw" & (ah? + ap? — 2vh) Fi (Pn a) 
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That (4.15) is truly a solution of our problem can be proved assuming that the 
function f(r) is such that the above series and the series obtained by applying the 
2 


operator Z and converge adequately. 


62” 

When v = 0 (i.e., 7 = 3) the problem (4.13) consists of finding the bounded 
steady temperatures wu (r, z) in the cylinder r € a,z > 0, if it is assumed that heat 
transfer into surroundings at temperature zero takes place through the surface r = a, 
according to the linear law u, (a, z) = — hu (a, z). 


Remark 4: The problem (4.13) is usually solved by means of the finite Hankel 
transform only in the case n = 3 (Colombo‘, p. 82). Now, by combining the finite 
transforms (4.1) and (4.6), it is feasible to solve this problem for each n => 3, even 
more, for an arbitrary integer 7. 
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L’-CONVERGENCE OF A MODIFIED COSINE SUM 
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We introduce here a new modified cosine sum and study its L'-convergence 
to a coine trigonometric series belonging to the class S of Sidon5 


1. INTRODUCTION 
Let 


co 


es =f a, cos kx AAD) 


k-1 


be a cosine series satisfying a, = 0(1),k + oo. If there exists a sequence <A,> 
such that 


Ax 1 0, k cee pan lee) 
S Ap < so 03) 
k=0 
Of Aue ARR (1.4) 


we say that (1.1) belongs to the class S introduced by Sidon’. 


Let the partial sum of (1.1) be denoted by S, (x) and f(x) = lim S, (x). 


Noo 
Concerning the L!-convergence of Rees-Stanojevi¢é cosine sums‘ 
3 Ses Ak k 
=}+2tAa a, cos kx 
Sa (x) + Fe k + ed Sion 


to a cosine trigonometric series, belonging to the class S, Ram! proved the following 
theorem : 


Theorem A—lf (1.1) belongs to the class S, then 
If — fall’ = 9 (1), 2 — ©9. 


In the present paper, we introduce a new modified cosine sum as 
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g(x) = 5 Ls > 


k=1 ja 


Ss A (4) k cos kx 
k 


and study its Z}-convergence. 
2. LEMMA 
The following lemma will be used in the proof of the theorem : 
Lemma 1*—If | cx, | < 1, then 


TT 


\|> cp EDA dx < Cont 1) 


0 k=0 2 sin 
where C is a positive absolute constant. 


3. RESULT 
We prove the following result. 


Theorem—Let (1.1) belongs to the class S. If lim | a,4; | log n = 0, then 





n—->oo 
lf — gnll = 0 (1), n > 0, 
PROOF: We have 
— 4% ak aks. a, 
(= 7+ 5 kos kx(a (F )+a Geena: (*)) 


k=1 


pee Qk Ana 
D + > boos kx | % rad 


k=1 
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= ee dS ax cos kx — ath Ss k cos kx 


| | k=1 


Gir TF, 
= S,(x) - creat D,, (x). 


Now, making use of Abel’s transformation and Lemma 1}, we have 


(ise if > Aax Dx (x) | dx 


0 0 k=n+1 


(equation continued on p. 1103) 
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0 
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Under the assumed hypothesis, x(k + 1) MA, converges and therefore the first term 
in (3.1) tends to zero asn > ©. Moreover, by Zygmund’s Theorem’ (p. 458). 
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since f | D, (x) | dx behaves like log n. 
=r 


The conclusion of the theorem now follows from (3.1) and (3.2). 


Corollary—If (1.1) belongs to the class S and lim | a,,, | logn = 0, then 
Noo 
IF — Sail = 0 (1), 2 > ce. 


PRoor: We notice that 


SIF) - S,@) lax = TIS) — 8.) + 8 &) — 5,0) | dx 


< f1f@) ~ 8) | de +18, (@)-S, (0) | dx 


TT 


< firey scolar | 





eee D: (x) dx. 





n+ ] 
-1 
wT wT 
Since lim | | f(x) — g, (x) | dx = 0 by our theorem and | Gni1_ py (x) dx 
PER n+] * 
-oF ef 


behaves like | a,,1 | log n by Zygmund’s Theorem cited above, for large values of n, 
the conclusion of the corollary follows. 
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The hydrodynamic stability of an annular liquid jet having a regular or 
irregular mantle solid axis, subjected to capillary and inertia forces, is present- 
ed. An eigenvalue relation valid for all modes of perturbation is derived, 
using the energy principle. The characteristics of the model are identified 
analytically, confirmed numerically and interpreted physically. The model 
is stable to all non-axisymmetric modes but it is unstable only to axisymmetric 
sausage) modes whose wavelength is longer than the circumference of the 
liquid jet. However the maximum temporal amplification values prevailing 
of such model are far lower than that of the full liquid jet. The greater the 
radius of the cylindrical solid axis; the slower the corresponding growth rate 
values but the greater oscillation frequency values in the stability regions i. e., 
the thicker the solid mantle, the larger its stabilizing effect. The present 
results reduce to those of Rayleigh'* if we impose that the radius of the cylin- 
drical solid axis tends to zero. 


1. - INTRODUCTION 


The experimental and theoretical hydrodynamic stability of a full liquid jet has 
been treated comprehensively since more than a century. This is not only from the 
academic viewpoint but also for its crucial applications in miscellaneous domains. 
Research in this field intensified when it became apparent that the physical properties 
of liquid jets play a fundamental role in a rapidly growing number of applications such 
as: the spinning of synthetic fibers, fuel atomization, spray drying, the production of 
controlled surfaces for heat and mass transfer in industrial and engineering processes 
and even the diagnosis of certain abnormalities of the human Urinary tract. 


It was Plateau’s observations® about the stability of a liquid jet (subjected to the 
curvature pressure) which led him to attribute the capillary instability to the surface 
tension. He was the first to obtain the critical wave length experimentally and theore- 
tically using a naive approach. The decisive break through came with Rayleigh’® who 
devised an elegant mathematical model for the breakup of liquids jets. Rayleigh” laid 
inexorable foundations for the theoretical treatment of such problems and developed 
the most important concept of the maximum mode of instability. By extending 


Rayleigh’s theory, Weber" considered the capillary instability of a viscous liquid jet. 
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These and other various problems are summarized by Rayleigh’® and, later on, also by 
Chandrasekhar’ but with miscellaneous extensions. 


The effect of non-linearities on the capillary instability of a hydrodynamical jet was ex- 
amined by Yuen'®, Wang", Nayfeh*, Nayfeh and Hassan’ and a complete analysis was 
finally given by Kakutani ef al.? using the derivative expansion approach developed by 
Kawahara‘. 


The capillary instability of an annular liquid jet (a liquid jet having a vacuum or 
gas Core jet) is also important to be investigated. The response of an incompressible- 
inviscid annular liquid jet subjecting to surface tension and the inertia forces (as the 
inertia force of the liquid is paramount over that of the gas jet) was given by Chan- 
drasekhar’ (p.540) for the axisymmetric mode only. Recently Kendall® have performed 
very interesting experiments with modern equipment for the capillaryinstability of an 
annular liquid jet for all modes of perturbations. Radwan!® investigated analytically the 
capillary instability of that model in incisive account, and it is found that the theore- 
tical results are in good agreement with Kendall’s experimental results’. Indeed, 
Kendall® explained clearly and speculatively, in a large context, about the important 
applications of that model. Moreover he attracted the attention for the theoretical 
studies of the annular liquid jet in general. 


The endeavours of the present work investigating the hydrodynamic stability of 
an annular liquid jet having a circular solid axis as a mantle subjecting to the capillary 
and inertia forces, by employing the energy principle. 


The results of the present work reduce to those of Rayleigh’ if we impose that 
the radius of the solid axis tends to zero. 


2, FORMULATION AND EIGENVALUE RELATION 


We consider an incompressible-inviscid liquid jet of radius R coaxial with a solid 
axis of radius R, (= qgR where 0 < gq < 1) in equilibrium state. A cylindrical coordi- 
nates (r, p z) system will be used with the z-axis coinciding with the axis of the coaxial 
(solid-liquid) cylinders. The influence of the vacuum medium surrounding the model 
is inevitable and neglected as long as the velocity of the perturbed liquid jet is not too 


large. The equilibrium density of the liquid P is assumed to be uniform. The gravita- 
tional effects are totally excluded. 


To carry out the present analysis by employing the energy principle; we should 
find out the total kinetic energy EF and the total potential energy Vin order to write 


down the Lagrangian equation of motion. It may be noted that the Lagrangian fun- 
ction L is constructed as 


L=E-—-v walt 


and the equation of motion is 
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Zi aL Lip 
dt sax) aye ...(2) 


where a* is the Lagrangian variable for the problem at hand and dot over a* denotes 
the time derivative. 


Let the equilibrium be disturbed then, in a cylindrical polar coordinates (r, Q, Z) system. 


According to the linear theory, the perturbed (liquid-vacuum) interface will be described 
at any time by 


r= R(1+ Q) vate) 
with 


Q= 2 2 am (k) expi(k z+ mo). ...(3a) 


Here Q is the elevation of the surface wave normalized with respect to R and measured 
from the equilibrium position. The coefficients am (kK) are functions of time and are 
much smaller than unity so that their higher orders can be neglected, based on the 
linear perturbation technique. Moreover am(k) are complex and the condition that the 
elevation is real being 


a_m (—k) = at (k) (4) 


where the asterisk* over am (kK) implies complex conjugate. m and k are the azimuthal 
and longitudinal wave-numbers respectively where m is integer while kK may have all 
continuous (real -) values. Now since the motion is irrotational and the liquid is invi- 
scid, the perturbed velocity vector u can be derived potentially viz. 


u= grad dQ. ...(5) 
Using equation (5) with the equation of continuity for incompressible fluid, the potential 
velocity ¢ satisfies Laplace’s equation 

ve = 0. ...(6) 


Note that Euler’s equation of motion has been used only implicitly in taking the irro- 
tational flow as presistent. In view of the 9 and z-dependence (cf. equation (3)); 
equation (6) reduces to an ordinary differential equation whose solution is given in 
terms of Bessel functions with argument of purely imaginary values. Under the present 
circumstances, the non-singular solution for ¢ must be 


$ (r; @, 2; t) = = E [b,, (k) Im (kr) + cm (k) Km (kr)] exp i (kz + m Q) 
mk 
best a) 


where bm (k), Cm (k) are unspecified functions of time and /m (kr), K,, (kr) are the 
modified Bessel functions of order m of the first and second kind respectively. 
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Using eqns. (6) and (7); the total kinetic energy E can be expressed (on using Gauss 
theorem) as a surface integral in the form 


E = 3p §¢ (grad $. dS) 


od ad 
=tese(r Sar 2S 47) dade ...(8) 


lg 
where the surface elements d S(= (r do dz, dr dz, r d¢ dr)) in the orthogonal curvilinear 
cylindrical coordinates (r, p, z) have been used. Since ¢(r, 9, z; t) is of first order; 
ean. (8) can be rewritten, up to second order, as 


fa HeR a fa ee 9) 





Substituting for ¢ from (7) into (9) and integrating with respect to 9; we get 


E = pa RJ dz & EI (bm (k) Im (KR) + em (k) Km (KR) [bm (1) 


It, (R) + c_,, (1) K‘, (R)] exp? i(kz + m 9) ...(10) 


where the prime on Bessel functions denotes the derivative with respect to r and where 
k and / are different dimensional longitudinal wave-numbers. 


The relation among the coefficients am (k), bm (k) and ¢m (k) can be determined by im- 
posing the boundary condition that the normal component of the velocity must be 
compatible with the deformed (liquid-vacuum) interface (3) at r = Rand that vanish 
atr = R,. This yields 


d ¢ 





ae = 0 atr = R, = gR B.f48) 
or 
= af atr= R : tL) 
from which 
b,, (k) = — cm(k) KK’,  (qx){Ii, (qx) ...(12a) 


cm (k) = x-* R® am (k) Ti, (qx) U0), (gx) Ki (x) — Ki (qx) I, @)? 
...(12b) 


where x (=kR) is the longitudinal non-dimensional wave-number, 


= Now, following Chandrasekhar? (p. 539), the potential energy V of a system 
arising from capillary forces is simply proportional to the total superficial area S; i. e. 
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V=TS ai) 3) 


where 7 is the surface tension coefficient. For the deformation given by (3), the total 
potential energy being 


ow 


v-T\l4+( 50 e+ re st Jt yt dz dp 


& 


=TR | a: aQ+ou+ad oy ( 52 )r+ ae( 32 pr de. 


... (14) 


t 
S eee 
™ 


By a resort to equation (3) and that the liquid is incompressible (the volume is coserv- 
ed); equation (14) gives 


Vad RT+oRT| dz a3 SD (1+ FR +m) am (b) a (k) 


m k/ 
x exp i (kz + mo). OATS) 


Therefore, for a single mode (since there is no interference) given by 


Q = am (k) exp i(kz + mg) + a_m (—k) exp (—i (kz + mg)) 
= a, (k) exp i(kz + mo) + a® (k) exp (~i (kz + mg9)) .»-(16) 
the total potential energy (per unit length in the z-direction) is given by 
fe eRYT 2 — m= x) a, (k) an, (KDI. ...(17) 


Similarly from eqns. (10) and (12) and for a mode given by (16); an expression for the 
total kinetic energy (per unit length in the z-direction) can written down. 


By constructing Lagrangian function [cf. eqn. (1)], equation (2) gives following 
equation of motion 


0 = a, (k) = 7 ( ~ mm — x) a, (kb) 


x ( ba ie et (qx) — id (qx) Ki, (x) ' 


UN Pe ee 


( In (x) Ki, (x) — Ty, (2) Kn (x) ) ...(18) 


Now, as usual for stability problems; based on the linear perturbation technique, 
we assume that the time dependence is in the form exp (c f); where « the growth rate is 
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of the perturbation and if it is imaginary o = iw then w/27 is the oscillation frequency. 
Henceforth equation (18) yields, at once, the eigenvalue relation 


oc = 


a ee ..(19) 
= 5 Rell m x*) F,, (x) ( 
where 
x| 1,0) Ky (@x) — 1,4) Ky (x) | 


F(x) = 
[ tm (2) Ke) ~ L,, G2) Km (XD | 


ooe(19a) 


3. DisCuUSSION AND CONCLUSION 


Equation (19) is the desired eiegnvalue relation of an annular liquid jet having a 
circular solid axis as a mantle, subjecting to the liquid intertia force and endowed with 
surface tension at the (liquid-vacuum) interface. By means of that relation the charact- 
eristics of the present model can be determined: one can identify the instability regions 
(in particular their critical wave-numbers, maximum growth rate values and the cor- 
responding wave-numbers) and those of stability as well. 


The eigenvalue relation (19) relates the growth rate « (or rather the oscillation frequency 
«) with the entity (7/p R°)-)/2 as a unit of time, azimuthal wavenumber m, longitu- 
dinal non-dimensional wavenumber x, Wronskian like expression [cf. eqn (19a)], the 
geometric factor q (the radius of solid axis normalized with respect to that of the liquid 
cylinder) and the physical cylindrical functions appropriate to the present model. 

As a limiting case as q tends to zero; equation (19) reduces to 


2 


ot ae (1 — mt — x4) (x 1, (2m (2). ...(20) 


This is the classical dispersion relation of a full liquid cylinder, as was derived by 
Rayleigh; for its discussions we may refer to Chandrasekhar? (p. 537). 


Other limiting case can be considered here as q + 1; the stability characteristics of 
such a case can be discussed as follows. Using a series development for the modified 
Bessel functions around the point gx, and neglecting terms of order (1 — q)’, 


the eigen- 
value relation (19) for m = 0 yields 
3? I, (x) Ky (x) — Fx Ry 
——.=5 x? (] — x2 aS ae Oe aN, 
TPR! * U— XD aha Gy aye 
By the use of the well-known (Wronskian) relation 
W (1m (X) Ky (X)) = Im (x) Ki, (&) — I, (x) Ky (x) = — x71 onl aa) 


equation (21) gives 
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TpR® = x° (1 — x) (1 — 4). ...(23) 


From equation (23) it is clear that o is depleted when g $ 1, so that the instability 
Sh hbaad slower and slower, the thicker the solid cylinder is with respect to the liquid 
jet; as is intuitively clear. Taking the derivative of o° we get 


d 
Fe 2 (TIP RY (V— gh (L — xt (1 — 22%). -.(24) 


From eqn. (24) we obtain directly that the maximum value of the growth rate at 
x, =1/(2)'!? = 0.707; that agrees perfectly with the value determined numerically and 
is barely more than the x,, (= 0.697) for the full liquid jet. 


By an appeal to the recurrence relations (cf. Abramowitz and Stegun’) of the 
modified Bessel functions 


2 1 (x) = Imus (X) + Iman (X) ...(25) 


2K’ (x) = — Ky-1 (x) — Kms (*) ...(26) 


and for each non-zero real value of x that Im (x) is always positive and monotonic in- 
creasing while K,, (x) is monotonic decreasing but never negative; we can observe that 


I’ (x) is always positive while K,, (x) is always negative. According to the foregoing 


arguments; one can show that 
EF (xy>.0 rat2t) 


for each non-zero real value of x, all q values and for all modes of perturbation; and 
that Fm (x) is never change sign. By a resort to the inequality (27); eqn. (19) yields that 
o? < O forall m+ O, but that o? > Ofor—1<x< land o2 <0 forx >1 or 
x < — lifm =O. Therefore, the model is stable to all non-axisymmetric modes, but 
is unstable to axisymmetric (sausage) modes whose wave-length \ = 2z/k is longer than 
the circumference 2 R of the liquid jet. Henceforth, for an annular liquid jet having 
a circular solid axis as a mantle we conclude that it is unstable only in the sausage 
mode in the domains x > 1 or x<—1; which is exactly the same for the full liquid jet 
subjected to the same forces as here. Indeed, this can be interpreted physically as 
follows. For such a model i.e. full liquid jet with or without solid axis; the potential 
energy due to surface tension (which is the only source of energy) can lower its value 
in the sausage mode m = 0 and hence give rise to kinetic energy, in particular there an 
enhancement in the very long wavelengths. That is also true even when the annular 
liquid jet having a vacuum or gas-core jet as a mantle (i.e. a hollow jet), see 
Chandrasekhar’ (p. 540). 


However for the present model, even if the solid axis is not regular we expect that 
the domains of instability are the same as if itis regular. One can make reasonable 
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estimates of the upper and lower limits of the growth rate as a function of x by consi- 
dering the thinnest and thickest radius of the solid, say 9g, R and q, R respectively. 
Then one may expect that the dispersion curve for the model will be in between those 
of cylinders with values.g, and q, respectively. Since the influence of g is not too large 
(see the numerical calculations Fig. 1), especially when g < 0.7 the upper and lower 
bounds will usually yield a good estimate if g, and g, do not differ too much. However 
if the internal solid has a periodic shape along its axis, there might turn up an enhanced 
effect in the dispersion relation for the corresponding wavelengths. 


To identify more clearly the effect of g-values on the instability characteristics of 
the present model, the eigenvalue relation (19) with m = 0 (since there is no instability 
in the modes m 0) and x > 0; has been used in the computer simulation for different 
values of q and for all (short and long) wave lengths. Thetnumerical results are present- 
ed graphically, see Figs | and 2 for instability domains respectively. There are many 





0.1 0.2 O34 2 06s S05 06 07 08 09 10 


Fic, 1. The eigenvalue relation for the capillary instability of an annular liquid jet having a 
circular mantle solid axis (in the sausage mode m = 0 with x <1). The abscissa measures 
the wavenumber in the unit 1/R and the ordinate the growth rate in the unit (T/eR°)}/2, 


features in these figures. Although the domain of instability (mainly 0 < x < 1 for 
m = 0) remains the same for all q-values; the area under the instability curves decreases 
with Increasing g and rather proportionally over the whole domain of instabilit 

Similarly we may prove for the stability curves (cf. figure 2), 
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YT/pR3 


10 
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Fic. 2. The eigenvalue relation for the capillary stability of an annular liquid jet having a circular 
mantle solid axis (in the sausage mode m = 0 with x 2 1). The abscissa measures the 
wavenumber in the unit 1/R and the ordinate the oscillation frequency in the unit 


(T/pR*)'?. 
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Moreover from these numerical results we deduce that the larger x is, the smaller 
is the decrease in sc and » with increase in g. However the effect is very small unless 
one considers very different wavenumbers. This result isin agreement with the very 
small shift in x,, from 0.697 to 0.707 for g increasing from 0 to unity. Particular at- 
tention is paid to the limiting cases g = 0 and g = 1. A somewhat large number of 
digits were needed to determine the wavenumber corresponding to the maximum 
growth rate. In fact with four relevant digits for ¢ it was still not clear whether or not 
the x-values corresponding to the omax was the same for all g. Using greater precision 
will clarify that x-value corresponding to omax does shift from 0.697 to 0.707 as q in- 
creases from zero to unity. Therefore, one can state that the solid backbone has a 
stabilizing influence on the liquid: its influence is more stabilizing in case of larger 
values or q, i. ¢. the larger the ratio of the solid axis to the liquid jet. However the 
stabilizing effect is rather small: for a large q like 0.9, the o-values decrease by a factor 
of about 2.3 and for a fairly extreme value like q = 0.99 they decrease by a factor of 
only 7.0. One can give a rough explanation for that weak stablizing influence as follows. 
As stated above the amount of energy sets free fromthe potential energy (the energy 
due to the curvature pressure) is the same for all g-values, but for a larger g (a relatively 
thicker solid cylinder), the amount of liquid requiring motion is smaller; thus the 
instability is not restrained very much in spite of the stabilizing effect of the solid which 
makes motion difficult. 


For values of q exceeding say 0.95, we may use the approximation formula (23). 
However, it is worthwhile to mention here that when the fluid layer becomes rather 
thin several effects occur which may influence seriously the present analysis. One effect 
is that the solid surface may be uncovered, yielding an additional amount of energy 
as mentioned above. Other effect is the Maragoni effect which is due to fluctuations in 
the surface tension and which is relatively more important for thin layers. 
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This paper contains the two-dimensional solutions of the problem in bipolar 
coordinates. The inclusions are elastic and perfectly bonded to an infinite 
elastic plate under the action of an isolated force. (i) ¥ applied at the origin 
in the positive direction of x-axis, and (ii) Y applied at the origin in the 
positive direction of y-axis. 


The stresses on the boundary, in each case, are found numerically and pre- 
sented graphically. It is found that the boundary «, = 0.8, the maximum 


values of the stresses 2x and 23 are (i) 15.10X units attained at 8 = 109.5° and 
10.47X units attained at 8 = 180° respectively, in the first case; and (ii) 
3.13Y units at 8 = 0° and 25.64Y units at 8 = 15° respectively, in the second 
case. 


INTRODUCTIONS 


The problem of stress distribution around elastic inclusions in an infinite plate 
has been studied by many authors. Saleme' has studied the stress distribution around a 
circular elastic inclusion in a semi infinite elastic plate under tension parallel to the 
straight edge. Mahata” has studied the stress distributions around two equal circular 
elastic inclusions in an infinite elastic plate under the action of (i) a centre of pressure 
radiating from a point, and (ii) a bi-axial or uni-axial tension. 


In the present paper we have studied, in the absence of the body forces, the stress 
distribution around two equal circular elastic inclusions in an infinite elastic plate 
under the action of an isolated force (i) X applied at the origin in the positive direc- 
tion of x-axis, in section 1, and (ii) Y applied at the origin in the position direction 
of y-axis, is secion 2% 


In each section, we find two solutions, one of which is regular inside and the other is 
is regular outside the inclusions. Numerical calculations are performed and graphs of 


stresses are drawn. 


The bipolar coordinates? are defined by 
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x+i(y+ a) aria 
e+ Bole FTG a Se? 
where a is a positive length. 
*. x = sin B/h, y = sinh a/h, ah = cosh « — cos P. .s(1,2) 


Now, « = constant, represents a set of co-axial circles having the two poles A (9, a) 
and B(o, — a) of the above transformation, for limitimg points. The curves 
8 = constant are a system of circles through A and B and interSecting the first net 
orthogonally. > 0 on R.H.S. of y-axis and 8 < 0 on its L.H.S. while on the y-axis, 
8 = 0, except on the segment AB where 8 = +7. At infinite, « = 0, 8B = O and at 
A, B, « — cc and « + — co respectively. 


The stress function %, in the absence of the body forces, satisfies the differential 
equation, 








ot a4 a4 e2 e2 oy, 
Camere wel meer Oe 
ae CL) 
The stresses in terms of AX are 
oe h ss inh lays apo yey 
aaa = [(cosh « — cos 8) ae — sinha — sin op 
sates) 
~~ o2 
aa = — (cosh « — cos 8) ba op (™)- sil 1S) 
The components of the displacements are 
7 ({ @ sinh « _9 sinh 
pape Sey ( éa cosh a — ae Nas ( @B cosh «—cos a) 
(AQ) ...(1.6) 
ey Pye) me bot sin B : ( es sinh « 
tains a »( 8B cosh « — cos 8 ) (Ax) + “6a cosh «—cos 3) 
(hQ) sila) 
where AQ is the associated displacement function and satisfies the differential equation 
3? oe o2 
én op (#2) = ( opie BAR 1 ) (hn. ...(1.8) 


METHOD oF SOLUTION 


Let the inclusions be defined by 


Oy aS ee ey ee 


Pa 14 & 
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Section I 


For an isolated force X applied at the origin in the positive direction of X-axis, 
X 

Xo = — ——— (8 — vx log r) A ip Bok 
7 


where (x, y) and (r, 9) are the cartesian and polar coordinates of any point in the 
plate. 


: in peas Ss lp. _1{ sinha) _ 7 
ht. = ee | sin “% tan ( aT) 4vsin Bloga 
L¢ cosh « + cos B 
tv sin B log (£28 cosh « — cos B )] 
X wT < : 
= | - 5 sinh « + Dye (x) sin n8 APSE) 
1 
where for « > 0 
] 
Ki) = —(5 + 1)e™+ 1+ vloga | 
and for tide ; : 
; S$ a v > 
Birer2, Kaige 62 {1° 1} ao mente ke | 
+ cosh «)]. J 
...(2.4) 


For the complete stress function we have to add to %, a stress function %; which 
will give no stress at infinity and satisfy eqn. (1.3). We assume, 


es as na ..(2.5) 
Ss Ore ) 


where 
¢, (%) = A, cosh 2« 


a ¢, («) = A, cosh (n + 1) a + B, cosh (n —1)a,n > 2 | ...(2.6) 


The associated displacement function 4Q, corresponding to AX, is given by 


co 


a, = « > + n (%) cos nf eetacay 
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where 
ti (4) = —2A, sinh 20 | .«-(2.8) 
v, (2) = —2 A, {sinh (n+1) « + B, sinh (1 —1) «,} nD2. 
The complete stress function is given by 

hk = hk + hh. re FN) 
The stresses ac, op and displacements u., up corresponding to AX are given by 


_~ co 


eas DY [(1 — n®) cosh a {, (4) + Kn («)} — sinh « {4% (a) 
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+Ki (v)} +4 (n — 1) ( — 2) {bea (@) + Kys (@)} 


+ $(n + 1) (n + 2){bn4; (&) + Kyi1(«)}] sin nB — ...(2.10) 
aeB an ¥ Ta) Fk! (@) 4 [= cosh a (oo (oh ER ee 


+ K), («)}] cos n8 
oe Sy | - n cosh « {f' (a) + Ki (a)} 
2 
+2(n— 1) {d,_. (a) + KL, (@)} + 4 t+ 1) {6,, (@) 


+ K‘,. (a)} Joos nd (2M) 


Eh = 
XPna =a+y[ Fa — cosh « cos 8) + >, Ls (a) sin np | 
1 


+ >, (>) (8, @) cosh a — 44, (@) — 3 44, (@) 


~ sinh « $, (@)} +n yy (a) cosh a — 4 (n — 2) dy, (a) 


— 2 (nN + 2) Yagi (a)] sin 18 mA Ae) 


STRESS DISTRIBUTION 1119 


a pat Ue: a| Me (a) — Ss sinh « sin B + »S M,, («) cos nb | 





— {(1 — v) di (@) +4 4) (@)} + > [1 — v) {n cosh « 


¢, (x) a 3 (n rae 2) Pn-1 («) Ts 3 (n 2 2) Pnti («)} 
+ cosha ¥) (x) — 4 4, (@) $4, («) — sinh a y, ()] 





cos n8 (2:19) 
where Z; («) = 2e-* (sinh 22—1)— 2 e-?™ sinh « + vsinh « loga 
for 
no 2, CL. (a) = e™ [{(— 1)* — 3} (= — sinh «) + (— 1)" v sinh 2a 
4 oS [{( — 1)" + 1} cosh « (nm cosh « + sinh 
«) + 2v sinh « (nm sinh « + cosh a) ] 
M, («) = v (log a — 4 — 2e cosh «) + (1 — e73*) 
M, («) = v (—cosh « loga + 2e* — 8 e-2% cosh «) 
— (= , sinh « + e™® sinh 2a ) 
Mz (a) = > + Ne (a) 
forn > 2, 


Ny (a) = 2+ (=1)" (1 — v) sin 2a e-* 


— [2v cosh « (n cosh « + sinh «) 


na 


. e~ 
+ ((—1)" + J} sinh « (n sinh « + cosh «)) ~ 54 


and for n > 3, N, («) = Mrz (%). 


Consideration, analogous to those introduced earlier, regarding the single valuedness of 
the displacement field, leads us to take the stress function within the inclussions, in 


the form, 
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oo 
ike i Rano gees SS $,(a)sin m6] (2.14) 
1 
where 
Ds > } 
¢, («) = Are fora > 0 | 
= A, e* fora < 0 : 
and | 
forn > 2,¢,< A, e ("tie 4+ Be e-(-1) for a > 0 | 
= A, e*= 4 B eln- for a < 0, J Ath) 
The associated displacement function corresponding to AX is 
oo 
hQ = > dn (%) cos nB ...(2.16) 
1 
where 
v, (a) = 2 4, (a) 
and for 
n> 2, b, (a) = 2¢, (a). A213) 


The stresses within the inclusions are 


op = es > [((1 — n*) cosh « @, (a) + $(n — 1) (n — 2) ris (a) 


+ 4(n + 1) (n + 2) 6,41 («) — sinh « ¢/, ()] sin nB »..(2.18) 


co 


ab= He, + > GO-DH,O+E@ + DE, 


1 


— ncosh « ¢, (%)} cos n]. »-.(2.19) 


The displacements #,. and ag within the inclusions are for « > 0, 


Ehig ve A Kai . = Wf j tae, 
XPna ~ ~ 40(1 —v ) sinh « sin B + 4 Gy (a) — (1 — v) ¢, (a) 





(equation continued on p. 1121) 
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a at iy) in cosh ud, (4) — ¥ (1 — 2) $.xa (2) 


= 3 (n a7 2) ra (a)} 
—coshat, («) +4 O_, (@) +4 O,,(@) + sinha ¥, (a)] 


x cos nf ...(2.20) 


Ehits ~ 3 . y $ 
Wea = (1 — v) Ay (1 — cosh « cos B) + > [C1 — v) {cosh « 4, («) 


1 





— sinh « ¢, («) — 4 ran (x) — 4 4/4, («)} 


—ncosh « o, (a) + (n — 2) dy-, (@) + 4 (0 + 2) Una 
(«)] sin np. 


BOUNDARY CONDITIONS 


In case of perfect bond between the plate and the inclusions, we must have on 


a=+ 4%, 


oct = an, a8 aa. pekarl) 
Ua. = te, uB = ug. Bol 4 


By eqns (2.10), (2.1 1), (2.18), (2.19) and (3.1), we have 
4 Uy (ay) = bn (oy) = bn (01) + Ky (a1) aGo) 
XE d, (4) = by (1) = (4) + Ky (4), ~..(3.4) 


Ci kl Pee ee 


By eqns. (2.12), (2.13), (2.20), (2.21) and (3.2) we have 


(1 + y) 7 = (i— y) Ay ...(3.5) 


Ay = [1 +9) Mo (as) + Ge (l= 9) Ks (es) — Ki IMS >) 
(3.6) 
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3A, V3 (4) + B, {3V, (%,) — a’ sinh a1} — A, {2V, (a) cosh «, 
E eT ‘ 
+ a’sinh a1} — Ful — v) {Kj (%) cosh a1 + $ Ky (a1) 


+ K, (%,) sinh «,} + 2K, (%,) cosh «, — 3K, («,)] — (1 + v) 
sae BN (a,) = ( al 


A, V; (a,) + Bs ig («,) + a’ cosh %,} —A, ae (%,) cosh «, + 2V; 
, E Xj 
x (%1) — a’ cosh a} + E [U1 — v) {K, (x1) cosh «, 


— 2K, (%1)} — 2K; (a) cosh a + Ki (01) + 2K; (a) 


sinh (%,] — (1 + v) M, (4) = 0 ...(3.8) 
and for n > 2, we have, 
(1 + 2) Vag2 (21) Ants + {(m + 2) V, (%,) — a’ sinh n a} Br+ 
— A, [2n Vay; (%1) cosh x, + a’ sinh nx] — B, [2n cosh x, V,_1 («,) 
— @ sinh nx] + A,_; [(n — 2) V, (01) + a’ sinh na] + (mn — 2) 


x Vacs (22) Bona + pr (1 — 9) (Ky (4) cosh 24 — 9K: (2%) 

eels Kis, («;) — sinh «, K, (%1)} — 2n cosh % K, (a1) + (n + 2) 

x Ky-1 (01) (n +2) Kags («,)] — (1+ v) LZ, (a1) = 0 ...(3.9) 
Vans (3) Ans) ot {Vi («;) + a’ cosh na} Basi = [2 {cosh ay ary (a1) 

— sinh % V,,; (%,)} — a’ cosh nx;] A, — [2 {cosh x, ie (x) 

— sinh | an” (x)} + a’ cosh na] B, or [V,, («,) — @ cosh na] 

X Ay + VO Veep ee eo v) 

n n—-2g \%1) Bn-1 E — v) {n cosh a K, (%1) 
= 4 (n = 2) Reet (%1) <= 4 (n — 2) Kyat (x,)} ar 2cosh aT Ki (%) 


as K,-1 («) + Kis, (4) +2K,, (%,) sinh %J—(1 + v) M, (%,) = 0 
ons Szh0) 
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where 
a’ =] -v— (1 ay) 


b=1+yv— (1-7) 


eu 
Vz (a) = — 5 sinh mz + cosh nx 


b 
Vv, (a) =n (— 7 cosh nx + + sinh nx), 


Thus eqn. (3.5) gives A,; eqn. (3.6) gives A,; eqns. (3.7) and (3.8) give A,, Bs; and eqns. 
(3.9), (3,10) give A,, 8, forn = 3, 4,...: Finally the equations (2.10) and (2.11) give 
the stresses outside while eqns. (2.12) and (2.19) give those inside the inclusions. 


Section II 
For an isolated force* Y applied at the origin in the positive direction of y-axis, 
gi 
% = ee (x8 + vy log r). .-- (4.1) 


Since the stresses and displacements corresponding to the stress function cosh x —cos 8 
are zero, sO any constant multiple of it is omitted from the stress function, and finally 
we havefor « > 0 


7 


co 
[7 Pe 5 [sinh « (1 + v log a) + 5 sin 8B + >. G, ©) c08 n6} 
1 





...(4.2) 
where 
G, (4) =v(1 — e-*) — 1 
and for | Fs eae | 
fee eV er ar | 
Js a [e-(@n-)% en (nt lay, | ».(4,3) 
n J 
The associated displacement function AQ, is given by 
oo 
AQo  _ 9 Ss G,, (x) sin n8. ..(4.4) 
Y/27 
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Here we assume, 
co 
ih, == ae [Bo ~ (cosh « — cos B) + DS ¢, («) cos nB] ...(4.5) 
7 
] 
where 
¢, («) = C, sinh 2a 
and for ; .. (4.6) 
n> 2, ¢,(%) = C, sinh (n + 1) « + D, sinh (n — 1) « 
co 
hole 5-128, a cosh wt > b, (@) sin 7B] AAT) 
1 
where 
1 (x) = 2C, cosh 2« 
and for 
n> 2, v, (%) = 2[C, cosh (n + 1) «+ D, cosh (n — 1) a]J° ...(4.8) 
The complete stress function is given by 
hy. a hXo “+ hx. ...(4.9) 


The stresses «x, «8 and the displacements uo, up corresponding to AX are given by 


Fie = — By sinh « (cosh « — cos 8) + ¢, («) + G, («) 
+ “S [(1 — n*) cosh « {9 (2) + G, (a)} + (2 — 1) (n — 2) 
x {$n-1 (%) + Gry («)} + 4 (n + 1) (nm + 2) {bay («) + Gasr (%)} 
— sinha {$/ («) + Gi («)}] cos nB .»- (4.10) 
eo — By cosh « sin B + ea) sin 28 +> [n cosh « {’ (x) 
1 


+ G, (@)} —3(n— 1) OF, ®O+G4_@}—ta@4) 


x {bin, () + Gi,, («)}] sin nB ates 
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GEhue 
Yj2n_ ~ Bolt (1 — v) — (1 + v) cosh? « + 2v cosh acos 8 + 4 (1—v) 





x cos 2B + 28 cosh « sin 8] — 3 (1 — v) $, (a) + ti (x) 


]+y 


+ 5 Gi (x) + (1 — v)(1 + v log a) (1 — cosh « cos 





B) — nae — v) sinh sin B + 5 ((1 — y) [cosh « 
1 


{p; («) + G* («)} — sinh « {¢n (%) + Gr (~)} + 4 {¢; (a) + G.,, (~} 


— 4 {9n4, () + Gi, (@}H + n cosh « {— 4, (a) + 2G, («)} 


F (+ 2) {8 Yaar (%) — Gras (@)} + (2 — 2) ha — Gn-, («)}) 
X cos n8 <oah4.22) 


aEhug : 
7 BoB sinh « cos 8 + (1 — vy) [— _ (1 + cosh «) 


7 
+ —5 cosh « cos 8 = ae (cosh « — 1) cos 28 — sinh « 


(1 + v log a) sin 8] 
+ 8 — v)[— ncosh « {¢, («) + Gy (x)} + $(n — 2) 


X {n-1 (&) + Gra (%)} + 4 (2 + 2) (bgt (%) + Grai(%)}) + cosh 


a {i (x) — 2 G («)} + sinh « {—v, (a) + 2G, («)} 
+{-34 #, @+ Ge NO) tel 3 Vi. () + G4_, (@)} 


sin np. ...(4.13) 


Within the inclusions, we take 


hi = pt | Be sinB + S ‘ («) cos nB] ] ...(4.14) 
1 
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where 


and for 


where 


and for 
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= } 
¢, (a) = C, e-™*, fora > 0 | 
= C, e* fora < 0 | 
b. | 
n> 2, 6, (a) = C, ete 4+ D, e(*-1*, fora > 0 | 
= C, elmtia 4 Dy elm" fora <0 J (4.15) 
. aie 2 Non ee 
* hO'= >. v, («) sin nB ...(4.16) 
1 
di (4) = — 24, (~) 
n>2, U,(«) = — 26, (a) 4. r=s(4.17) 


The stresses and displacements corresponding to hi are 


a, 


eaee é, (~) + SS: [1 — n*) cosh « g, (%) + 4 (mn — 1) (nm — 2) bp-y (2) 


+ 4(n + 1) (n + 2) ony, («) — sinh ag’ («)] cocnB. — ...(4.18) 


~ 


i 3 (n cosh x $, (2) — 4 (n— 6, (®)—4FN +1) S,, 





(~)] sin nf. ...(4.19) 


ahE iin l+y-, > B S 
viet =a t (4) + (1 —v)[- B, sinh «sin? + > {cosh « 
1 





x $,, (a) — $9, , (2) — 44, (2) ~ sinha ¢, (a)} cos nB] 
+ z [— cosh « Un (ct) + $(n — 2) Vn-a (#) 
+  (n — 2) daar («)] cos m8 ...(4.20) 


ahEiig 7% D . 
yaw ~ (1 — v)[— Bo (1 — cosh « cos 8) + 3 {—n cosh « $, (2) 





+ 41 — 1) bys (@) + 4 (n + 2) bays (@)} sin 8] 


(equation continued on p. 1127) 
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< > [cosh ah) (a) -— 4 V, (~) —3 U,, @) — sinh « 


Y, (2)] sin n8- (4.21) 
Boundary Conditions 
We have on the boundary « = + «,, 
Ux = Ue, Up = ip. tS 2) 
Therefore from (4.10), (4.11), (4.18), 4.19) and (5.1) we have, 
By ae 0 (Sea) 
bn (41) = bq (41) + Gp (a) 
ae! 
4’, (as) = 6, (2) + Gi (mr) Cy 
for te 233. .; 
From eqns. (4.12), (4.13), (4.20), (4.21) and (5.2), and using (5.4), we have 
=o -— E 
By = (1 — vy) a/2 (1 — v) = ... (5.5) 
1 — v) (1 + vlog a) 
C= —- [2 Gi (a1) + fa CF vig |/cos at, #456) 
l+v— Pp (1 + v) 


3C, (Cy — a’) cosh 3a, + Dy (3Co — a’) cosh a, + 2C, (a’—C, cosh 2%) 
x cosh «, + (2a’ — Cy) cosh a, G, (#1) — 2 sinh a G, (a) 


— 3Cy Gp (a1) — a’ Gp (a,) — 2 (1 — v) (1 + v loga) cosh «, = 0 
5.7) 


3C, (a’ — C,) sinh 3a; + D, (3a’ — Cy) sinh «, + 2C; [(Co — a’) cosh 


4 ) bas Co ) 
a, sinh 2a, + Cy sinh «,] + 3G, («,) (a 4 


1 


—34(1+ s) G5 (#1) — 2a’ cosh «, G, («,;) — Co sinh a G' 


(equation continued on p. 1128) 
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(a;) + Co cosh a, Gi (a,) —2(1 — v) (1 + v log a) sinh «, =0 
.. (5.8) 
and for n > 2, we have, 


; n+2 ,, : 
bce. (a) Caan + Des aT capes Ri (%;) + a’ cosn na, | 





ae [= fa cosh a, Rss (a,) + a’ cosh na ] 


nl 


are tS [== cosh ~, R’_, («,) + a’ cosh nas | 


n—2 





+ Cy4 Ri, («;) — a’ cosh na, | “Diy Ron (hel 


+ a’ | cosh a, Gi (a,) — sinh «,. G, (a,) — 4 Gi_, (@,) 


= 4 G54, (4) | + Colm cosh 2 G, (a) - 24? 6, (md) 


—_ = Gai (a) | = 0 ...(5.9) 


(n + 2) Rage (41) Cog + Dagi [2R, (%1) — a’ sinh noi] — C, [a’ sinh Na, 
2 sinh «, R’ 





+2 (n+ 1) cosh «, R,,1 (a) — Fra n+} (a:)] 
ears ate 2 sinh « - 
n [2 (a — 1) cosh # Ry-1 (a1) —-—>——*- Ri_, (a) 


— a’ sinh nz,] + C,_, [nR, («,) + a’ sinh no] 


n—2 
yi 





+ (n — 2) R,-2 (%,) D,-1 + a’ [n cosh «, G, («,) — 


n+2 
x Coy (a1) — ae pose Co (a,)] + Gs [cosh Ay tps («,) 


— sinh 2, G,(%)— 4G,_, (a) — 46%, (u))}=0 ...(5,10) 


where C, = 2(1 — E/E) 


R, (a) = 3 sinh no + ¢ cosh nx 


; b 
Ri (a) =n lees cosh na + sinh na), 
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Thus (5.3) gives Bo, the equation (5.5) gives By, the equation (5.6) gives C,, the equa- 
tions (5.7), (5.8) give C,, D.; and equations (5.9), (5.10) give C,, D, for n = 3, 4, ... 
Finally, the equations (4.10), (4.11) give the stresses outside and eqns. (4.18) and (4.19) 
give those inside the inclusions. 


NUMERICAL CALCULATION 


Numerical Calculations are carried out with a = 1,v = v = }, EJE 3, «, = 0.8. 
The circumferential stresses are found Table | and are shown graphically for both the 
cases against 6 (Fig. 1, 2). 


STRESS 


STRESS 











Fic. 1. Stresses on the boundary «,; = 0.8 Fic. 2. Stresses on the boundary «; = 0.8 
under the force X along the x-axis. under the force Y along the y-axis. 
TABLE I 
Max. value For isolated force X parallel For isolated force Y-para- 
of the stress to x-axis llel to y-axis 
aa 15.10X at 8 = 109.5° 3.13Y at B = 0° 
ap 10.47X at ® = 180° 25.64Y at 8 = 15° 


a ESE 
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The effect of periodic variation of suction velocity on free convection flow 
and heat transfer through a porous medium is in vestigated. The problem 
becomes three dimensional due to variation of suction velocity in transverse 
direction on the wall. A series expansion method is used to get the solution 
of the governing equations and the expressions for velocity and temperature 
fields are obtained. The skin friction and the rate of heat transfer at the wall 
are analysed in detail. 


INTRODUCTION 


The problem of laminar flow control has gained considerable importance in the 
field of Aeronautical Engineering in view of its application to reduce drag and hence 
the vehicle power requirement by a substantial amount. The development on this 
subject has been compiled by Lachmann’. Theoretical and experimental investigations 
have shown that the transition from laminar to turbulent flow which causes the drag 
coefficient to increase, may be prevented by the suction of fluid and heat transfer from 
boundary layer to the wall. Gersten and Gross? have studied the effect of transverse 
sinusoidal suction velocity on flow and heat transfer along an infinite vertical porous 
wall. The flow and heat transfer through a porous medium, however, has not received 
much attention despite its application in many branches of engineering. In this note, 
the effect of suction velocity variation on free convective flow and heat transfer in 
a porous medium is investigated. 


ANALYSIS 


We consider free convection flow of a viscous incompressible fluid through a 
porous medium bounded by an infinite vertical porous wall. A coordinate system 


with the wall lying on x — z plane and y-axis perpendicular to it and directed into the 
fluid is introduced. The suction velocity distribution is taken in the form : 
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Vy (Zz) = v, [1 + € cos | --(1) 


which consists of a basic Steady distribution VY, < 0 with a superimposed weak 
transversally varying distribution « V cos nz/L of wave length ZL. Since the wall is 
infinite in x-direction, hence all physical quantities will be independent of x, however, 
the flow remains three dimensional due to variation of suction velocity. Let us denote 


velocity components u, v, w in X,Y, Z directions respectively and temperature by T. 
The governing equations are 








a 4 = 0 Wel 
re tw Ho BIT —T.) Hy (Se + ar)- eA) 
y+ rape tS )-% (4) 
v 5 th at ~ «(SS + ae) bat) 


In eqn. (3), variation in density is taken into account only in the derivation of 
the buoyancy force while other density variations are neglected within the frame-work 
of constant property fluid. The last terms in eqns. (3), (4) and (5) account for the 
pressure drop across the porous material. It is worthwhile to mention that the basic 
flow in the medium is entirely due to buoyancy force, caused by temperature difference 
between the wall and the medium. 


The boundary conditions of the problem are : 


y=0:u=0,v=y,(z), w=0,T=T, 


are 
yoo:u=-0,w=0,p=p.n,T=T x. (7) 


When the amplitude ¢ of oscillations in suction velocity is small, we assume 
the main flow velocity u in the following form 
u=umt+eu, teu,+... ...(8) 


: T — T 
The similar equations hold for other variables v, w, p and @ = ey ae : When 


e = 0, the problem reduces to two-dimensional free convection flow in an infinite 
porous medium with constant suction velocity at the wall. In this case eqns. (2) to (6) 


reduce to 
% _ 9 at) 
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d dfuy _ We 

% Be = gf (Tw — Too) 80 ss dy? K 
d0 wt d* 6, 

ray a Rapes 


The solution of these equations is 


Gy 


ee ka 


(e~™Y —_ e- PRY), Vo = Vo, Wo = 0 


Bo = ePRY, Dy = Doo 


where 


--.(10) 


oh ann 


(iz) 


iehke) 


2 aD 3 
7 = 'y|L, P = va, Goede Me tcc pe pee es 


y2 





LV 
My fo ee 


Rie 


iF ae 
are ay) 


+ (R*4 + K,) 2, 


When «€ + 0, the series expansion ($) is substituted in equations (2) to (6) and 
like powers of € are equated to get the perturbation equations of various order in «. 
For small values of ¢, it is sufficient to consider the perturbation equations only of 


o (€), which are 


ov, ow, 











oy ar Patten 





20, 26, (86, , 8%, 
mtn Beno 4 Bey 


with boundary conditions 


y= 0:u,=0,», =Vy cost 7, = 0,6, = 0 


y > co:u, = 0, w, = 0, p; = 0, 6, = 0. 


| 


..(14) 


wi I5) 


...(16) 


acta) 


...(18) 


(19) 


This is a set of linear partial differential equations which describe the three-dimen- 


sional cross flow. 
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First of all we shall consider the eqns. (14), (16) and (17). The solutions for 


vi (¥, Z), W1 (¥, z) and p, (y, z) are independent of main flow component u, and 
temperature field 6,. We now assume v,, w, and p; as 


Vs, (9, 2) = Vo m v1, (9) cos nz ...(20) 
Wi (9, Z) = — voviy (}) sin xz ites) 
Pi(}, 2) = PVE pis (9) cos xz ten) 


where z = z/L and prime denotes differentiation with respect to J. 


Equations (20) and (21) have been chosen so that the continuity equation (14) is 
satisfied. Substituting these expressions into (16) and(17) and applying the corres- 
ponding transformed boundary conditions, we get the values of v,, w, and p, as: 








(9,2) = os d (xe — de-™) cos 72 et eed 

Wi (J, 2) = —t (e 47 — e- 7?) sin nz (24) 
/ A = 

Pi (9, 2) = Nz PVG Ee=ed e-™” cos 7Z ate) 


where 


K, eas 2 2 1/2 
N,=1+ ph andA = 5 (RYU4 + n* Kt, 


Assuming u; and 6, of the form 
U, (9, Z) = uy, (J) cos mz ..(26) 
61 (J, Z) = O11 (9) Cos mz lai) 


and substituting in eqns. (15) and (18), we get 


y weG TR silo ERP A As ) e-¥7 
a cee (arn ~ any — At Art 4s Je 





— r = - a 
ss -(Atm)y _ _“_ g-(rt+m)y _ 4, e~(AtPR)Y 
TKN, 7Ns 
+ A, e-(w+PR)Y + Ag em | cos 7Z .--(28) 


e~4y ] COs 7Z 


tee) 


. P* R Tens born) aeteRD As 
6; (9, Z) ar aoe Ee e-(AtPR)y ap? P3 


where 
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P, P? 
RA (P=) eek 


~~ 
Il 


A(1 + P) + KR, Ps = 


N,; = Lk (R? + 4K,)/3 
m 








a (1 + PP,/RP,) 
Ms = 1+ K,/2ma, A, = cee oa 
Acme os AP IR A Pi 
*" P? R?+ 2a PR — RP — aR — K, 
7 PR 2 2 R2/4yij2 
As = P3(n/P, — AlnP), A = a E(w? + PAR?/4)P. 


RESULTS AND Discussions 


We now discuss the important flow characteristics of the problem. The expres- 
sion for shear stress along the direction of x is obtained as 


Gy 

aa Tx = i \ dy a) 

Cee Pv/L? p v?/L2 

= - (PR — m) + «GF, (P, R) cos xz ...(30) 
1 
where 
R 7m Ane 
F, (P, R) = Bie a | aN, + Mee + m— A) 


+ A, PR + A, (A — PR — 7) + A3(A -%)] mo es 


The skin friction factor F, (P, R) in the limiting cases, R>Oand R > eco, tends to 
zero. Equation (31) is numerically evaluated for different values of R and permeabi- 
lity parameter K; and the results are shown in Fig. 1. From Fig. 1, it is obvious that 
F, tends to its limiting values as R > 0 and R > co for each K;. It is also seen that 


F, incre h it decreases and 
as the value of K, is increased, 


ases with R until it attains a maximum value, after whic 
approaches to zero. It has also been observed that, 
the skin friction factor F, decreases significantly. 


The heat flux at the wall in terms of Nusselt number Nu is given by 
—~ dw k (72 ) 
Nue= ———w Kf OB 
PVy cp (Ty — To ) PVocp \ay JyLo 


= 1+ e[! — F,(P, R)) cos 22 ...(32) 


where 
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Fic. 1. Variation of skin friction factor F, with Fic. 2. Variation of heat transfer factor F, with 





Rfor different values of permeability Rfor different values of permeability 
parameter K,. parameter Kj. 
] A nP APR 
F, (P, R) = Aj —-— — —- —-_ — 
, (P, R) &iPp(G as = (A + PR) : “|. 
..(33) 


It is found that F, takes constant values for limiting values of R. When R > 0, 
F, > 1 which is obvious because there is no oscillatory fiow. However, when R > ©9, 
F, > 0 which shows that heat transfer approaches to quasi-steady value. These 
results are shown in Fig. 2. It is interesting to note that F, (P, R) increases with 
permeability parameter K,. 
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MHD SWIRLING JET WHICH ORIGINATES FROM 
A CIRCULAR SLIT 
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(Received 21 November 1986; after revision 27 June 1988) 


The effects of an axial magnetic field, which varies inversely as the radius 
vector, on the velocity distributions in a swirling jet of a viscous incompressi- 
ble electrically conducting fluid are studied. It is noted that the radial as well 
as tangential velocity decrease near the slit of the jet with the increase in the 
value of the magnetic interaction parameter, These decelerated fluid particles 
move in the positive axial direction and the points where they exactly balance 
the motion of the incoming fluid have been calculated for different values of 
the magnetic interaction parameter. 


NOMENCLATURE 


B = magnetic field vector 
Bo = axial magnetic field at a unit radial distance 


Fy function introduced in eqn. (3.2) 





Gj = function introduced in egn. (3.3) 
Jo = linear momentum flux at a large distance from the slit 
Z, = initial angular momentum flux 
Se B; : 
nm = p (magnetic parameter) 
Q = volume flux in the radial direction 
r = radial coordinate 
u = radial velocity component 
y = axial velocity component 


w = circumferential velocity component 


Z = axial coordinate 
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Greek Letters 
( 3J, 1/3 
a = ————$— 5 
4m Pr/v 


> In \iP 
p= (3) 





— = non-dimensional variable ( ee ) 
rV/yv 
= coefficient of dynamic viscosity 
= density 
vy = kinematic viscosity 
os. = electrical conductivity 
vb = _ Stokes’ stream function. 


1. INTRODUCTION 


The flow of a laminar jet issuing froma slit (plane free jet) and a jet issuing from 
a circular orifice (circular free jet) without swirl, of an electrically non-conducting, in- 
compressible viscous fluid were investigated by Schlichting!, Bickley® and the corres- 
ponding MHD flow has been studied by a number of workers including Peskin’, Smith 
and Cambel!, Pozzi and Bianchini®, Bansal’, Gupta’ and Mishra and Bansal®. 


In swirling flows, besides axial and radial velocity components, the tangential 
component of velocity should also be considered. It is both the linear momentum and 
the angular momentum of developing swirling flow which play an important role in 
determining its ultimate form. The case of a laminar swirling jet of a non-conducting, 


viscous, incompressible fluid was considered by Loitsianski’, Gortler!®, Steiger and 
Bloom" and solutions were obtained using the assumption of similarity of the velocity 
profiles. Similarity in turbulent swirling wakes and jets have been studied by 
Reynolds'® and Chervinsky’’. Chervinsky noted that for the turbulent axisymmetrical 
swirling jets, the asymptotic solution is valid only in regions where the axial pressure 
gradient is negligibly small. 


Recently, Mishra and Bansal’* studied the effect of circular magnetic field on the 
decay of a weak swirl in the axially-symmetrical circular free jet of an electrically con- 
ducting, viscous fluid and found that the swirl velocity increases due to the presence 
of the magnetic field, though its decay is inversely proportional to the distance square 
along the axis of jet, as is the case in the non-magnetic flow. 


In the present paper, we have studied the effects of an axial magnetic field, which 
varies inversely as the radius vector, on the velocity distributions in a swirling jet of a 
viscous incompressible electrically conducting fluid which originates from a circular 
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slit. A perturbation on the Loitsianski model is applied and first order perturbation 
solutions for the velocity distributions are obtained. It is observed that the radial as 
well as tangential velocity decrease near the slit of the jet with the increase in the value 
of the magnetic-interaction parameter. These decelerated fluid particles, because of 
the law of conservation of mass, move in the positive axial direction and the points 
where they exactly balance the motion of incoming fluid have been calculated for 
different values of the magnetic-interaction parameter. 


2. FORMULATION OF THE PROBLEM 


Let a rotating incompresible, viscous, electrically conducting fluid be discharged 
through a circular slit formed by two circular discs of negligible radii and negligible 


distance apart in the presence of variable axial magnetic field B, (0, 0, Bo r-') ‘and!mix 
with the same surrounding fluid being initially at rest (Fig. 1.). Taking the origin in 






TY 
Se SY 


Bisis: cia 
r M 


w 


Fic. 1. Spread of MHD swirling jet which originates from a circular slit. 


the slit, the boundary layer equations in cylindrical polar coordinates governing the 
fluid motion may be written as : 


du du w? 07x mu 


uy 


or 02: aie One en Sal stack) 
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ow ow uw 0°w mw 
u :- + y TPS + — v aor = Saree »»-(2.2) 
é ; 7 
be (ru) 4- see (rv) = ee 
where 
ce BS | 
/(i—— p isi(2.4) 


The boundary conditions are 








Ou ow 
Ba horas me Osi Opa 0 
z>boru>—,w> = ; AcE 


where compatibility conditions are used in obtaining the conditions at infinity. 


Besides these boundary conditions the following integral conditions, which govern 
the linear momentum flux (at a large radial distance from the slit) and the initial angu- 
lar momentum flux of the jet respectively, should also be satisfied (Loitsianski'*). 


co 
lim 27Pr | u. dz= J, tte. 0} 
r—+>0oo ties, 
Qn p a Oia Wo dz a Lo ater) 
—0o 


where up, vo and w, are the velocity distributions when the magnetic field is zero 
(m = 0) and are given by 











are = SAL ABS 
eee ta tS op Or (2.8) 
by avy Eri fu (2.9) 
Wp = Bat z pet vis (E) ...(2.10) 
PR ear aah 
ce nv/v 
18 _)” (2.12) 
‘ ( 4m Pv 
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= (5) SIPAE). 
fe = tanh (E/2), f,, = 0. 
£11 = 0, 812 = 8/2 sech® (&/2). ...(2.14) 


3. ANALYSIS 


The equation of continuity (2.3) is identically satisfied if we consider the Stoke’s 
stream function U,,, such that 


BPs ert ee ena waaay 
Oz or 


The momentum equations (2.1) and (2.2) can be reduced to a set of ordinary differen- 
tial equations if we consider the following series expansions for /,, and w, which are 
perturbations on Loitsianski’s model’®; i. e. equations (2.9) and (2.10) 


Yn = % fv cS (gist 1A (Ai os S S a = ) Fy «© | A) 
s=6 o=) 71 
and 
v=se[> replay + 3 Ss aA Gj © | mes, Re Fat 


f=] 


where Fj; and G,; are unknown functions of & to be determined. 


Using (3.!) the expressions for u, v and w from (3.2) and (3.3) are obtained as 


See ee im? wy 
u=* Ge to he hide (F;, + —F he Jab | ] 








...(3.4) 
Aha/ F 2 ‘ 
y= vl ( ef — fio) + rt iene 
+ DER, — Fo + FL ++ “4 Re) 
and 
a l 
ie a (811 + p S12 +e - (Gy, Ts - Gi, + gas] ...(3.6) 


where a prime denotes differentiation with respect to &. 


Substituting (3.4) to (3.6) in equations (2.1) and (2.2) and equating the coefficients 


of power of (m/a*) and then of r-?,, we get the following set of ordinary differential 
equations for firstord er perturbation : 
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Bivcthe Fs, 12h, +E’, fie Fi. ...G.7) 
Btenerat sf, Fk. = —2(¢,, 6. 4G, Gu), ...(3.8) 
Se hG. =O. 3.9) 
eget Cs ot fg Cig =i815 (1 — FY) — &1, Fro. ...(3.10) 


The corresponding boundary conditions from (2.5) are 


be Oc Fy = 0, Fi, = 0; F,, = 0, Fi, =.0; Gi 0,G! 203,11) 


E=co: F), = 1: Fi, =0;G,, = 0;Gi = B. s«s( 3-212) 


From the boundary conditions it follows that the solutions of equations (3.8) and (3.9) 
are 


F,, = 0, G,, = 0 er era hc 
whereas the solution of eqn (3.10) may be obtained as 


Giz (5) = B Fi, (&) (3.14) 


Hence, we have to integrate only eqn. (3.7), under the boundary conditions (3.11) and 
(3.12). It may be noted that it is an exact differential equation and integrating it we 
get 
g 
(1 + cosh &) Fi) = C (€ + sinh &) + | (1 + cosh &) log (1 + cosh &) d& 
0 
ato. Lo) 


where the constant C is obtained by the boundary condition F,, (eo) = 1. 


The integral in the R. H. S. of (3.15) can not be obtained easily in a compact 
form and therefore should be evaluated numerically for different values of &. Instead 
of this we have preferred to integrate the inhomogeneous third order ordinary linear 
differential equation (3.7) numerically by standard techniques using Runge-Kutta Gill 
method. It is found that 


Fy, (0) = — 2.52406 SAAT, 


and the results are tabulated in Table |. 
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TABLE I 
ee i eee 
E fio tae fe Fio F is 
0.0 0.00000 0.50000 0.00000 0.00000 —2.52406 0.00000 
0.2 0.09967 0.49503 —0.04934 —0.50080 —2.46416 0.59318 
0.4 0.19738 0.48052 —0.09484 —0.97310 —2.29127 1.11962 
0.6 0.29131 0.45737 —0.13330 —1.41102 —2.02433 1.52667 
0.8 0.37995 0.42782 —0.16255 —1,78337 — 1.69056 1.78524 
1.0 0.46212 0.39322 —0.18172 —2.08483 —1.32039 1.89209 
i ee 0.63515 0.29829 —0.18946 —2.51253 —0.41160 1.64611 
2.0 0.76159 0.20999 —0.15993 —2.53508 0.27420 1.08509 
7 Mes 0.84828 0.14021 —0.11894 —2.28474 0.68553 0.58700 
3.0 0.90515 0.09035 —0.08178 —1.88374 0.89189 0.26806 
35 0.94138 0.05691 —0.05357 —) 41257 0.87891 0.10024 
3.8 0.95624 0.04281 —0.04093 1211529 1.00030 0.04738 
4.0 0.96403 0.03532 —0.03405 — — — 
4.5 6.97803 0.02173 — 0.02125 — — _ 
5.0 0.98661 0.01329 —0.01312 — — — 
J 0.99186 0.00811 —0.00804 — — — 
6.0 0.99506 0.00493 —0.00491 — — a 





LS 


Finally, if we confine ourselves to the first two terms in the series expansions, we 
find 








u poe cg Wiss 
(ose = 2 Yes (&) + a? 10 (E)] ohSt?) 
yV . m ee 
(a +/v/r) Chis — fio) + ou? (EFI) — Fro) (3.18) 
and 
; ee @® + Gs Ol=- oe 
(Wo)r—0 8 81» wi 12 ( )] ms COE ...(3.19) 
The volume flux in the radial direction is given by 
eo 
gee | mae = 4nar /vT fio (eo) + = Fi) (ce)] 
a = m 
4nar Jv [1 — 1.11529 oak ...(3.20) 
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Fic. 2. Radial and transverse velocity distribution of a MHD swirling jet which originates 


from a circular slit. 


This shows that the volume flux in the radial direction decreases linearly with respect 


to ma]*. 
in the axial direction will be exactly balanced 


The points where the incoming flow 
yv=0 


by outward flow of the decelerated fluid particles is obtained from (3.18) when 


i. €. 


+0-2 


io 
4 5 6 
m - 
oe 
0:15 
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Fic. 3. Axial veolcity distribution of a MHD swirling jet which originates from a circular slit. 
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Sw (&) — Ef? (€) 


m 
2 





R 


OF ay hie (&) (3.21) 


4. NUMERICAL DISCUSSION 


As the magnetic field can only retard and not produce flow in the reverse direc- 
tion, the velocity field and the volume flux should be non-negative. Therefore the range 
of the magnetic-interaction parameter inthe present series solution should be 0 < 
m]a* <_ 0.198, a range in which the fluid has always a forward flow. 


The non-dimensional velocity profiles for radial and transverse velocity compon- 
ents, for various values of the magnetic interaction parameter (m/«”), are plotted 
against the similarity variable in Fig. 2 and for the axial velocity in Fig. 3 respectively. 
The effect the magnetic field is to decrease the radial as well as tangential velocity near 
the slit of the jet and to increase it far from the slit. From Fig. 3, it is clear that the 
axial velocity increases algebraically throughout the velocity field with the increase in 
magnetic-interaction parameter. 


Following physical explanation may be given for the results obtained : 


Due to the action of the centrifugal forces the fluid near the slit of the jet will be 
thrown outward and to compensate this a flow in axial direction towards the slit will 
follow. Now, the Lorentz force retards the radial and tangential motion of the fluid 
and therefore these decelerated fluid particles, from the law of conservation of mass, 
start moving in the positive axial direction and thus affect the incoming axial velocity, 
which moy eventually change its direction. It may happen that for a prescribed value 
of m/a? the incoming flow is exactly balanced by the upward fiow. Such points may be 
obtained from Fig. 3, where v is zero. 
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Eratta 


“Corrections to The neighbourhood number of a graph” 
by P. P. KALE AND N. V. Desupanpe, Indian Journal of pure appl. 
Math. 19 (1988), 927-929 


It needs correction as the result ny) < 2n, is false. The last proposition and its 
froof should read as follows. 
Proposition : For any graph G, 
Ny eae oy nerk0) 
Proor : Let T = {x,, X., ..., Xx} be a line dominating set of minimum cardina- 


k 
lityandx,;=mv,l <ick. LetH = U {w,v}. As every line in G is incident 
t=] 


with a point in H, H is a neighbourhood set of G, Hencemy < | H| <2|T7| = 2y’. 
This proves the proposition. 
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